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00 1 ABSTRACT 

o . 

. We explore the behavior of the blue galaxy fraction over the redshift range 0.75 ^ 

' z ^ 1.3 in the DEEP2 Survey, both for field galaxies and for galaxies in groups. The 

primary aim is to determine the role that groups play in driving the evolution of galaxy 
~ Q colour at high z. In pursuing this aim, it is essential to define a galaxy sample that 

I 1 does not suffer from redshift-dependent selection effects in colour-magnitude space. 

We develop four such samples for this study: at all redshifts considered, each one is 
complete in colour-magnitude space, and the selection also accounts for evolution in the 
galaxy luminosity function. These samples will also be useful for future evolutionary 
studies in DEEP2. The colour segregation observed between local group and field 
samples is already in place at z ~ 1: DEEP2 groups have a significantly lower blue 
fraction than the field. At fixed z, there is also a correlation between blue fraction 
and galaxy magnitude, such that brighter galaxies are more likely to be red, both 
in groups and in the field. In addition, there is a negative correlation between blue 
fraction and group richness. In terms of evolution, the blue fraction in groups and the 
field remains roughly constant from z = 0.75 to z ~ 1, but beyond this redshift the 
blue fraction in groups rises rapidly with z, and the group and field blue fractions 
become indistinguishable at z ~ 1.3. Careful tests indicate that this effect does not 
arise from known systematic or selection effects. To further ensure the robustness of 
this result, we build on previous mock DEEP2 catalogues to develop mock catalogues 
that reproduce the colour-overdensity relation observed in DEEP2 and use these to test 
our methods. The convergence between the group and field blue fractions at z ~ 1.3 
implies that DEEP2 galaxy groups only became efficient at quenching star formation 
at z ~ 2; this result is broadly consistent with other recent observations and with 
current models of galaxy evolution and hierarchical structure growth. 
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1 INTRODUCTION 

One of the most striking characteristics of the galaxy pop- 
ulation is the well-known environmental segregation of the 
two main galaxy types: red, early-type galaxies with little 



ongoing star formation preferentially occur in galaxy groups 
and clusters, while blue, late-type galaxies with significant 
star-formation activity avoid such systems and preferentially 
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populate the "field'Q- This observation has been recognized 
as a key to understandin g galaxy formation an d evolution 
for more than fifty years l|Spitzer fc Baa dc 1951). 

There is now overwhelming evidence that the galaxy 
population in clusters has evol ved significantly wi t h red - 
shift down to the present day. iButcher fc Oemlerl (ll984T ) 
were the first to present evidence that the fraction of blue 
galaxies, /{,, in clusters increases strongly with increas- 
ing z — the so-called Butcher-Oemler (BO) effect. This ba- 
sic result — an increased incidence of star-forming galax- 
ies in distant clusters — has been replicated in numerous 
later studies using a variety of star-formation indicators , 
including cluster blue fractions (iRakos fc SchombertJ 1 1995 1; 
Margoniner fe de Carvalhol 120001; iKodama fc BoweiT l200ll ; 



Ellingson et all l200ll ; IMargoniner et ail 120011 ), 
line galaxy fractio ns dPoggianti et ail l2006t)~nd morpho- 
logical fractions rtOemler et ail 1 19971: ICouch et ail 1998; 
van Dokkum et al] |2000| ; iFasano et ail l2000l : iLubin et ail 



20021 ; iGoto et al.M2003l ). Such studies have also been ex- 



tended to less massive galaxy groups ([Allington-Smith et ail 
Il993l ; IWilman et aill2005l ; IMartmez et al.ll2006l ). which also 
appear to have proportionally more star-forming galaxies 
back in time. To be sure, there remain strong reasons to 
question the veracity of the BO effect as it was origi nally 
presented for galaxies in the cores of rich clu sters ( e.q.. lKool 
ll98Sl ; ISmail etal.lll998l ; lAndreon et al.ll2006l ). But it is now 
indisputable that clusters, on the whole, had proportion- 
ally more blue, star-forming, and morphologically late-type 
galaxies — in short, more star formation — at z ~ 0.5 than 
they do at present. 

It is tempting to conclude that this evolution is respon- 
sible for the substantial growth that has be en observed in the 
number density of red galaxies since z ~ 1 (|Bell et al.ll2004bl ; 
IWillmer et ailfeoOrj ; iFaber et aill200^ . But it is also impor- 
tant to note that, in addition to evolution within the group 
and cluster environment, the build-up of the red galaxy pop- 
ulation could also be driven by an increasing number density 
of groups and clusters from the growth of structure. We shall 
try in this paper to shed light on the relative importance in 
DEEP2 of these two channels for red-galaxy formation. 

A blue fraction that declines with time in groups and 
clusters is a natural consequence of h i erarchical schemes 



for galaxy formation dKauffmann 



Diafcrio ct al. 2001; iBenson et ail 



.19951 ; iBaugh et ailll99rj ; 

200 ll ) if one assumes that 



groups and clusters play a role in quenching star forma- 
tion. Since these systems form at relatively late times, their 
member galaxies at intermediate redshift will have had less 
time, on average, to "feel" the effects of group membership; 
thus, a smaller fraction of them will have ceased forming 
stars. Indeed, it has been shown (Ko dama fc Bowerl l200ll ) 
that the evolution of in clusters out to 2 ~ 0.4 can be en- 
tirely accounted for by considering the effects of hierarchical 
structure formation and a declining universal star-formation 
rate. For exactness in what follows, then, we shall use the 
phrase "Butcher-Oemler effect" to refer just to this effect — 
a decline in the typical blue fraction in groups and clusters 
that results simply from the increasing age of these systems 



Throughout this paper, we shall use the term field to refer to 
those galaxies that are not in groups (i.e., isolated galaxies), not 
the full galaxy population 



with time, and not from any evolution in the quenching 
efficiency of groups and clusters. Although this working def- 
ini tion differs somewhat from the effect originally reported 
bv lButcher fc Oemlerl (| 19841 ). it reflects usage that has be- 
come common in recent literature. It is also worth noting 
that, under this definition, the BO effect may be stronger 
in some types of systems than in others (in clusters than in 
groups, for example), depending on the accretion rates and 
quenching efficiencies involved. 

It remains unclear, however, which specific physical 
mechanisms are responsible for quenching star formation, 
over what timescales they act, and even whether these mech- 
anisms are peculiar to groups and clusters. A variety of 
mechani sms has been proposed, i ncluding ( 1 ) galaxy merg- 
ers re.o. lToomre fc Toomrell 19721). w hich o ccur primarily in 
galaxy groups (e.g., Cavalie re et al.lll992l ) and can trigger 
AGN feedback that sweeps merger remnants of their gas 
jSpringel et ail 120051 ); (2 ) close, high- veloc ity galaxy en- 
counters ( "harrassment" ; iMoore et ail fl99^ ) . which occur 
primarily in clusters; (3) ram-pressur e stripping of galax y 
gas by the hot intracluster medium (|Gunn fc Gotdll972l ); 
and (4) the less violent process known as "strangulation" 
in which gas accretion onto galaxy di scs is cut off, either by 
stripping of galaxies' gaseous haloes (|Larson et al.|[l98fj|) o r 
by feedback from low-luminosity AGN (jCroton et al.l l2006). 
and star formation ends when the remaining supply of gas 
has been exhausted. There is evidence that each of these 
processes is at work on some level; it may be that no sin- 
gle mechanism bears primary responsibility for galaxy evo- 
lution. Disentangling how and to what degree each of the 
above processes shapes the galaxy population requires de- 
tailed observations of galaxies over a wide range of redshifts 
and environments. 

The observed evolution is rather complicated in its de- 
tails, however. For example, the well-know n morphology - 
density relation f or galaxies in clusters (e . q., lDresslerlll980l ) 
and groups (e.g. iPostman fc Gellerl Il984l) ev olves strongly 
with increasing redshift ( Dressier et ail 19971 ). and there is 
evidence that the evolution in intermediate-density envi- 
ronments has occured more rece ntly than in high-density 
environments l|Smith et al.ll2005l ). It also appears that the 
morp hological mix of cluster early-types has changed with 
time l|Fasano et al ] |200fj| ; |Postman et aill2005l ). Complicat- 
ing things further, cluster blue fractions d isplay a substantial 
system-to-system s catte r at all epochs (Butcher fc Oemlerl 
19841; ISmail et al l 1 19981; IMargoniner fc de Carvalhol |2000| ; 



Goto et ail 120031 ; |Pe Propris et ail \2004) . This scatter 



can be ascribed in part to the observed trends be- 
tween blue fraction and cluste r mass, i.e., velocity 
dispersion (INewberry et ail 1 19881). inferred virial mass 
(IMartmez et ail 12003) , richness IMargoniner et ail l200ll ; 
IGoto et ail 20031; iTovmassian et ail 120 05). or total opti- 
cal luminosity l|Weinmann et al.ll2006al ); the correlation be- 
tween and mass a lso appears to evolve with redshift 
(Pogg ianti et aill2006l ). The scatter may also arise in part 
from a correlation betwe en fb and a cluster's degree of dy- 
namical relaxation (e.g.. iMetevier et aill2000l ). The scatter 
in ft might mask or enhance the redshift ev olution, depend- 
ing on the cluster selection c riteria us ed (INewberry et ail 
1 19881 ; lAndreon fc Ettiorlll999D . Indeed, ISmail et ail (1998) 
observe no BO effect out to z ~ 0.25 in a sample of luminous 
X-ray clusters (but see also iFairlev et aill2002l ). and many 
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authors have commented on the presence of individual mas- 
sive clusters at hi gh redshift with low fb and well- formed 
red s equences (e.o lKoolll98ll : lrIomeier et al.ll2005l ; lAndreonl 
2006). It is clear that a robust study of galaxy evolution in 
groups and clusters requires a sample that contains a broad 
range of systems and is uniformly selected at all redshifts. 

Until relatively recently, studies of galaxy evolution in 
clusters were limited to cluster samples pre-selected from 
optical imaging or X-ray maps; such samples are prone 
to selection effects that bias the galaxy sample as a func- 
tion of z. The advent of large, densely sampled galaxy red- 
shift surveys like the 2-degree Field Galaxy Redshift Sur- 
vey (2dFGRS) and the Sloan Digital Sky Survey (SDSS), 
however, has allowed selection of large, low-redshift clus- 
ter samples directly from the galaxy distribution in red- 
shift space, where redshift-dependent selection effects can 
be well understood and accounted for. Several authors 
have studied the pro perties of the clust e r galaxy popula- 
tion in the 2dFG RS (TBalogh et al l 12004 |Pe Propris et all 
20041) and SPSS jGotol l2005al lbl; IWeinmann et all l2006al 



Quintero et al.l [20061 ) . Within SDSS, ther e have even been 



detections of an evo lving in clusters (|Goto et alj 120031 ; 
Mart inez et al . 2006) and o f an evolving morphology-density 
relation (|Goto et alj|2004l). With the high- redshift DEEP2 
Galaxy Redshift Survey (|Davis et al.ll2004l . Faber et al. in 
prep.) nearing completion, it is now possible to extend these 
studies to z ~ 1. 

DEEP2 has now yielded a catalogue of several thousand 
groups and small cluste rs over a wide range of masses at 
z ~ 1 dGerke et al.ll2005h . A study of the DEEP2 sample by 
ICooper et alj (|2006ph has already revealed that, without ex- 
plicitly considering groups, the correlations between galaxy 
properties and the density of the nearby galaxy distribution 
are q ualitatively similar a t z ~ 1 to what is observed locally 
(e.g., iBalogh et all |2004| ; iHogg et alj I2004T ). although the 
relations differ in detail. Interestingly, recent studies have 
shown that these c orrelations evolve with redshift, becoming 
stronger over time (iNuiiten et aLll2005l ; ICucciati et aLllioOa 
ICooper et al.ll2006al ). and that the color-density relation in 
the local Universe can be ascribed al most entirely to mech - 
anisms acting in groups and clusters (jBlanton et al.| [2006). 
The primary goal of this paper, then, is to establish the role 
that groups and clusters play, if any, in driving the evolu- 
tion of the DEEP2 galaxy population. In particular, we shall 
explore the possibility that group and cluster environments 
are responsible for the buil d-up of red galaxies observed over 
the DEEP2 redshift ran ge |Bell et al.ll2004bl ; IWillmer et al.l 
l200d : iFaber et al.ll2006l ); by so doing, we will try to shed 
light on the physical mechanisms driving the evolution. Al- 
though this study will focus on the evolutionary effects of 
galaxies' membership in groups and clusters (that is to say, 
of their dark matter halo masses), there is important com- 
plementary information to be gained by studying the effects 
of the local density of galaxies. Such a study is undertake n 
in a companion paper to this one bv lCooper et al.l (|2006af ) . 
For now, however, we will measure fb for galaxies in groups 
and clusters, as a function of redshift, and compare it to 
that found in the general field (i.e., in the galaxy popula- 
tion outside of groups). We will also investigate the relation 
between f b and cluster properties in order to gain insight 
into the processes at work and also to understand and con- 
trol any systematic selection effects. 



The thoughtful reader may find it perverse to use the 
blue fraction to study the evolution of red galaxies, but fb 
has significant historical precedent in the study of cluster 
galaxy populations, so we use it for consistency with earlier 
studies. In any case, the possibly more natural red fraction 
statistic is simply f r = 1 — fa. We have chosen to consider 
galaxy colour rather than other indicators of galaxy type like 
morphology or [O II] line emission primarily because it can 
be measured accurately for the largest number of DEEP2 
galaxies, allowing for robust statistics. The DEEP2 ground- 
based imaging lacks the resolution necessary for morphologi- 
cal classification of most galaxies at z ~ 1, and there is HST 
imaging for only a small fraction of the sample. Emission 
line strength can be measured accurately only in spectra 
with sufficient signal-to-noise ratio; not all DEEP2 galaxies 
meet this criterion. Also, it has recently been shown that [O 
II] emission may not give a cl ear indication o f star forma- 
tion, especially in red galaxies (|Yan et al.ll2006l ). Regardless, 
there is a strong and relatively tight correlati on between [O 
II] equivalent width and galaxy colour(e.<?., IWeinerl [20051 ; 
ICooper et al1l2006bh . so the two statistics should give simi- 
lar results. 

We shall proceed as follows. In SJH we describe the 
DEEP2 survey and the DEEP2 group catalogue, and in Sj3] 
we discuss the construction of mock DEEP2 galaxy cata- 
logues appropriate for this study. We describe our galaxy 
selection criteria and incompleteness corrections and present 
the precise definition of fb used here in 331 Measurements 
of fb in groups and the field are presented in and we 
perform tests with the mock catalogues in § [6] We discuss 
the implications of our results for galaxy evolution models 
in SJ7] and in Sj8]we summarize our results and conclusions. 
Readers uninterested in the details of our selection meth- 
ods and robustness checks should skip to § § 14. f I [S] and [7] 
Throughout the paper we assume a flat, ACDM cosmology 
with Q,m — 0.3. 



2 THE DEEP2 GALAXY REDSHIFT SURVEY 
2.1 Details of the Survey 

The DEEP2 Galaxy Redshift Survey is the first large, highly 
accurate spectroscopic survey of galaxies at redshifts around 
unity. As of this writing, the main survey observations are 
nearly complete (> 95%), with spectra obtained for 49, 220 
galaxies in four fields using the DEIMOS spectrograph on 
the Keck II telescope. The survey covers a total of ~ 3 deg 2 
on the sky to limiting magnitude Rab ~ 24.1. The bulk 
of these observations are of galaxies in the redshift range 
0-7 z < 1.4. Full details of the survey will appear in 
the upcoming paper by Faber et al. (in preparation), but 
summarize the necessary information for this study is sum- 
marized below. 

The four fields surveyed were chosen to lie in zones of 
low G alactic extinction using the dust maps of lSchlegel et al.l 
(1998). Three-band (BRI) photometry was obtained for 
each field using the CFH12K came ra on the Canada- F rance- 
Hawaii Telescope, as described bv lCoil et ail (|2004al ). Each 
field is covered by several contiguous photometric pointings, 
which are a convenient way to group and intercompare re- 
sults. Three of the DEEP2 fields each cover an area 90' or 
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130' x 30' on the sky with two or three contiguous pointings. 
In these fields, galaxies are selected for spectroscopy using 
a simple cut in BRI colour-colour space that has been op- 
timized to select galaxies at redshifts z > 0.75. This cut 
efficiently focuses the survey on high-redshift galaxies: it re- 
duces the portion of the spectroscopic sample at redshifts 
z < 0.75 to roughly 10% while discarding only ~ 3% of ob- 
jects at z > 0.75 (Faber et al. in preparation). Within each 
CFHT pointing, galaxies are selected for spectroscopic ob- 
servation if it is possible to place them on one of the ~ 40 
DEIMOS slit masks covering that pointing. Slit masks are 
tiled in an overlapping chevron pattern using an adaptive 
algorithm to increase the coverage in dense regions on the 
sky, giving nearly every galaxy two chances to be placed on 
a mask. Fur t her d etails of the observing scheme are given in 
iDavis et al] (|2004 ); overall DEEP2 targets ~ 60% of galax- 
ies that meet its selection criteria. 

The fourth field of the survey, the Extended Groth Strip 
(hereafter EGS), covers 120' x 16' on the sky. A concerted 
effort is underway by a large consor tium of observing teams 
(the AEGIS team; [5 avis et al"1l2006l ) to obtain a wide array 
of observations of this field from X-ray to radio wavelengths. 
Therefore, to maximize the evolutionary information that 
will be available, galaxies in the EGS have been targeted 
for spectroscopy regardless of estimated redshift, and at a 
significantly higher sampling rate than in other fields: each 
galaxy in EGS has four chances to be selected for observa- 
tion rather than two. This means that galaxies selected for 
spectroscopy in the EGS constitute a superset of galaxies 
that would have been selected using the criteria of the other 
three DEEP2 fields; hence it is possible to create a high- 
redshift subsample of the EGS whose selection (including 
sampling rate) is identical to that of the rest of the DEEP2 
survey. We will include this subsample in the present study. 

At this writing, spectroscopic observations have been 
completed for all three of the high-redshift DEEP2 fields 
and for more than three-quarters of the EGS field. All 
DEEP2 spectra have been reduced using an automated data- 
reduction pipeline (Cooper et al. in preparation), and red- 
shift identifications are all confirmed visually. Rest-frame 
U — B colours and absolute B magnitudes are computed 
us ing the K-correction algorithm described in Appendix A 
of IWillmer" et al.l (|2006l ). This study uses data from all of 
the DEEP2 CFHT pointings for which spectroscopic obser- 
vations have been completed. In each pointing, the fraction 
of the spectra that yield a successful redshift is greater than 
70%. 



2.2 The Group Catalogue 

The details of the DEEP 2 gro up-finding procedure are fully 
discussed in iGerke et all ((200 5|); we have now applied this 
procedure to the full current DEEP2 dataset, which is sig- 
nificantly larger than that used in the previous paper. We 
summarize the salient points of the group-finding algorithm 
here. 

Groups of galaxies in the DEEP2 sample are iden- 
tified using the Voronoi-Delaunay Method (VDM) group 
finder , which was originally implemented bv lMarinoni et all 
(2002). This group finder searches adaptively for groups 
(bound, virialized associations of two or more observed 
galaxies) using information about local density derived from 



the Voronoi partition and Delau nay complex of a given 
three-dimensional galaxy sample. I Gerke et alj (|2005h cali- 
brated the VDM group finder using the mock catalogues of 
lYan et all (|2003l l. achieving the primary goal of accurately 
reconstructing the bivariate distribution n(o, z) of groups as 
a function of redshift and velocity dispersion for dispersions 
a ^ 350 km s" 1 . 

For the purposes of this paper, however, we will be con- 
sidering the properties of galaxies, rather than group prop- 
erties, so we will focus on somewhat different measures of 
success here. In particular, this paper studies properties of 
the population of galaxies within groups (the group sample) 
and of the population of isolated galaxies (the field sam- 
ple), so it is crucial to determine the success of the VDM at 
identifying each of these populations. B y testing the VDM 
group finder on mock DEEP2 catalogues, iGerke et al l ([20051 ) 
showed that the fraction of real group members that are suc- 
cessfully identified as such (the galaxy success rate, Sgai) is 
0.79. Conversely, the fraction of galaxies in the reconstructed 
group population that are actually misclassified field galax- 
ies (the interloper fraction, fi) is 0.46. In addition, 82% 
of field galaxies are correctly identified, while only 6% of 
the reconstructed field sample is made up of misclassified 
group members. Both samples (group and field) are thus 
dominated by correctly classified galaxies, but each sample 
is contaminated by galaxies from the other. Therefore, any 
differences between the group and field galaxies should be 
somewhat stronger in reality than what the VDM recon- 
structs. As will be discussed below in § 14.11 the contamina- 
tion of the group sample is particularly bad for groups with 
velocity dispersion below 100 kms -1 , so we will reclassify 
galaxies in these groups as field galaxies. 

Finally, it is important to note here that groups in the 
DEEP2 surve y are typi c ally o f modest mass. The mock 
catalogues of lYan et all (|2003h indicate that the bulk of 
the DEEP2 groups should have virial masses in the range 
5 x 10 12 < Af v ir < 5 x 10 13 Af© (200 < o v < 400 kms -1 ), 
with very few groups having M v i r > 10 14 Mq; this is con- 
sistent with the estimates of the minimum group mass de- 
rived from the a utocorrelation function of DEEP2 groups 
|Coil et al.l l2006f) . In what follows, then, any conclusions 
drawn about the properties of groups should not be taken 
to apply to rich clusters. 



3 MOCK CATALOGUES 

The study of groups and clusters of galaxies is fraught with 
unavoidable sources of systematic er ror. For example, it has 
been shown (|Szapudi fc S zalav 199^) that an error-free clus- 
ter catalogue is not achievable for an incompletely sampled 
galaxy distribution, even if the physical distances to the 
galaxies are known. Moreover, in a redshift survey the pecu- 
liar velocities of galaxies induce distortions in the redshift- 
space distribution of galaxies, mixing the positions of galax- 
ies in groups with the positions of isolated galaxies and mak- 
ing accurate group detection even more difficult. Because of 
these unavoidable sources of systematic error, robust con- 
clusions require that we test our methods on realistic mock 

gal axy catalogue s. 

IGerke et all (|2005h tested and calibrated t heir group- 
finding methods with the mock catalogues of lYan et all 
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These catalogues were produced by populating dark- 
matter-only N-body simulations wit h galaxies followin g a 
halo model prescription that follows lYang et al.l (|2003h . In 
particular, the mocks are populated using a conditional lu- 
minosity function, <f)(L\M), that assigns galaxies to a halo 
of mass M according to a luminosity function whose pa- 
rameters depend on M. The parameters of the model were 
chosen to be consistent with the two-poin t correlation func - 
tion £(r) observed in early DEEP2 data (ICoil et al.ll2004ri) 
and w ith the local from Peacock et al. |Peacock et all 
2001). These mocks were sufficient for the purpose of test- 
ing our overall success at reconstructing groups, but they do 
not provide information about galaxy properties aside from 
luminosity. 

In the current work we shall need to test our success at 
reproducing the colours of the group galaxy population; this 
requires mock catalogues that include that information. The 
mocks must also reproduce any dependences of the galaxy 
colour distribution on local overdensity, since that is the 
trend we aim to probe in the data. To this end, we assign 
colours to mock galaxies by drawing colours from galaxies in 
similar environments within the actual DEEP2 data. This 
method is similar to, but less sophisticated than, the one 
used by Wechsler et al. (in preparation) to create the mock 
SDSS catal ogues that were us ed to test the C4 group-finding 
algorithm i|Miller et alj|2005h . 

The procedure is as follows. Local galaxy density in the 
DEEP2 data is measured by computing the distance (pro- 
jected on the sky) to each galaxy's third-ne arest neighbor 
within 1000 km s _1 in redshift, as detailed in lCooper et al.l 
(2005). DEEP2 galaxies contaminated by edge effects are 
discarded, and the data sample is limited to galaxies in the 
range 0.8 < z < 1.0 to minimize the effect of the apparent 
magnitude limit. In the mock catalogue, local density is es- 
timated using the projected distance to the seventh-nearest 
neighbor (within a catalogue of an enhanced spatial extent, 
to avoid edge effects in the final catalogue). The seventh- 
nearest neighbor in the complete, volume-limited mock cat- 
alogue has been shown to be a reasonable analogue to the 
third-nearest neighbor measu red in the more spa rsely sam- 
pled magnitude- limited data (|Cooper et al.ll2005T ). We then 
divide both the mock and data samples into quintiles of lo- 
cal density and identify each quintile in the data with the 
quintile in the mock that has the same local-density rank. 
The colour distributions within each density bin in the data 
can then be used to assign colours to the mock galaxies. For 
example, the 20% of mock galaxies with the highest local 
densities will be assigned colours randomly drawn from the 
observed colour distribution of the highest-density 20% of 
DEEP2 galaxies. 

Before doing this, however, we also divide the sam- 
ples by luminosity, since DEEP2 galaxy col our is observed 
to be correlated with luminosity (see, e.g., IWillmer et al.l 
2006). For the DEEP2 data, we sort galaxies into four 
bins in absolute B band magnitude. To ensure that sim- 
ilar parts of the luminosity function are being considered 
at all z, we shift these bins with redshift to account for 
evolu tion in the typical galaxy magnitude M* . (|Faber et alj 
2006) found that this value evolves as M* oc Qz, with 
Q — —1.37; we apply the same linear function of z to our 
luminosity binning of the data. The bins include only galax- 
ies brighter than Mb — 51og/i + Q(z — 1) = —20, since 



the sample is incomplete for fainter magnitudes (see Fig- 
ure [2]), and each bin is 0.5 magnitudes wide, except for 
the brightest bin, which includes all galaxies brighter than 
Mb — 5 log h + Q(z — 1) = —21.5. The mock catalogue is 
divided into the same four bins, except that in this case the 
bins evolve with redshift as M oc Q m0 ckZ, with Q mQC k = — 1, 
the M * evolution parameter that was assumed in Yan et al.l 
(2003). Also the faintest bin in the mocks includes all galax- 
ies fainter than Mb + Qmock(z — 1) = —20.5, which means 
that all mock galaxies fainter than this limit will have the 
same colour distribution as DEEP2 galaxies in the range 
-20.5 < M B - Q(z - 1) < -20. This is the best that can be 
done using this procedure, since the DEEP2 sample is in- 
complete for fainter objects; in any event, most of the mock 
galaxies in this regime will fall below the DEEP2 apparent 
magnitude limit, so there is little practical effect. 

Having divided the samples thus into bins of local over- 
density and absolute magnitude, we then add colours to the 
mock population by considering galaxies in each bin sepa- 
rately. For each mock galaxy in a given luminosity-density 
bin, a real galaxy is selected at random from the corre- 
sponding bin in the DEEP2 data sample, and we assign 
that galaxy's rest-frame U — B colour to the mock galaxy in 
question. This procedure produces a distribution in Mb vs. 
U — B colour-magnitude space that matches the distribu- 
tion observed in DEEP2 reasonably well (see Figure [I}. Be- 
cause the procedure only uses DEEP2 galaxies in the range 
0.8 < z < 1.0 and because it does not divide the samples 
by redshift, any intrinsic evolution in the colours of DEEP2 
galaxies will not be reproduced in the mock catalogues. This 
is desirable, however, since it will allow us to confirm that 
our selection and group-finding procedures have not intro- 
duced any spurious evolutionary trends. 

Using the colours assigned to the mock galaxies, it 
is possible to the K-correction procedure described in 
IWillmer et ail (|2006l ) to assign each galaxy an i?-band ap- 
parent magnitude, which can be used to select mock galax- 
ies with the same R < 24.1 apparent magnitude limit that 
is used in DEEP2. In addition to this magnitude limit, we 
also apply the DEEP2 DEIMOS slitmask-making algorithm 
to the mock catalogue, as projected on the sky, removing 
those galaxies that would not be targeted for observation. 
Finally, we dilute the remaining galaxy sample to reproduce 
the ~ 70% redshift success rate of the survey; the dilution 
procedure accounts for the slight magnitude dependence of 
this rate. 



4 SAMPLE DEFINITION AND 
MEASUREMENT METHODS 

4.1 Defining the galaxy sample 

Because the DEEP2 survey extends over a broad redshift 
range (Az ~ 0.7), selecting galaxies according to a simple 
apparent magnitude limit introduces selection effects that 
depend strongly on a galaxy's rest-frame colour. For exam- 
ple, the DEEP2 Rab ^ 24.1 magnitude limit corresponds 
to selection in the rest-frame B band at z ~ 0.7 and the 
rest-frame U band at z ~ 1.2. Therefore, galaxies with in- 
trinsically red colours will fall beyond the selection cut at 
lower redshifts than those with intrinsically blue colours. 
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Figure 1. Colour-magnitude and colour-redshift diagrams for the 
DEEP2 survey and for the mock catalogues described in the text. 
To aid comparison, the mock sample has been randomly diluted 
to contain the same number of galaxies as the data set, and each 
of the two data plots has the same contour levels as the corre- 
sponding mock plots (contours are spaced evenly in density). The 
subtle correlation between colour and luminosity in the data is 
reproduced well in the mocks at Mg — 51og/i < —20; at fainter 
magnitudes it is absent by construction. 



This effect is readily apparent in Figure [2] which shows 
rest-frame colour-magnitude diagrams (Mb vs. (U — B)) 
for DEEP2 galaxies, divided into redshift bins of width 
Az = 0.05. In each panel, distinct red and blue galaxy pop- 
ulations are apparent, with loci that are roughly divided by 
the dotted lines. The sharp cutoff in the galaxy population 
on the right side of each panel is caused by the DEEP2 ap- 
parent magnitude limit; this selection cut becomes increas- 
ingly biased against red galaxies as redshift increases and the 
observed R band moves further blueward of the rest-frame 
B band. 

It is obviously necessary to take this effect into ac- 
count when studying the redshift evolution of the blue frac- 
tion. In particular, we must ensure that galaxies of a given 
colour have been equally well sampled at all redshifts be- 
ing considered. The simplest selection method is to pro- 
duce a volume-limited catalogue with an absolute magni- 
tude limit — i.e., a vertical selection cut in colour-magnitude 
space. For DEEP2, however, this method severely restricts 
either the redshift range probed or the number of galaxies 
selected. For example, as can be seen in Figure [2] a limit- 
ing absolute magnitude of Mb ~ 51og/i = — 20 would allow 
measurement of ft only out to redshift z ~ 0.9; beyond this, 
the red galaxy sample would be incomplete, resulting in a 
spurious sharp rise in fb at all higher redshifts. On the other 
hand, a limiting magnitude of Mb — 5 log ft = —21.7 would 
give a complete, volume-limited sample out to z = 1.3, but 
such a sample would contain far too few galaxies for a robust 
measurement of /(,. 

However, to examine the evolution of galaxy properties, 
all that is required is to select a region of colour-magnitude 
space that is uniformly sampled by the survey at all red- 
shifts of interest. Such a selection cut is shown by the dashed 
curves in each panel of Figure [2] which are described by the 



(1) 

+ mm {[a(U - B) + b\, [c(U - B) + d]} , 

where zn m (equal to 1.3 in the figure) is the limiting red- 
shift beyond which the selected sample becomes incomplete, 
a, b, c, and d are constants that depend on zn m and are de- 
termined by inspection of the colour-magnitude diagrams, 
and Q is a constant that allows for linear redshift evolution 
of the typical galaxy lum inosity L* , as already mentioned 
in § [3] iFaber et alj l|2006f ) have measured the evolution of 
the galaxy luminosity function using data from COMBO- 
17 (|Bell et all l2004bf ) and DEEP2 at z ~ 1 in conjunc- 



tion with \aw-z data from the 2dFGRS (Madgwic k et al.l 
| 2002l; iNorberg et all |2002l) and SDSS l|Bell et al l 120031 ; 
blanton et all 120031 ') ; they find that changes in L* are well 
described by a linear evolution model, L* oc Qz, with 
Q = —1.37. By including this evolution in our selection cut, 
we are selecting a similar population of galaxies with respect 
to L* at all redshifts. 

We use Equation [1] to define three different samples, 
called samples II, III and IV, which are complete in colour- 
magnitude space to z = 1.0, 1.15 and 1.3, respectively; the 
parameters defining these samples are given in Table [1] In 
addition, we create a sample, called sample I, that is purely 
absolute- magnitude- limited (relative to M*) and complete 
to z = 1.0 by applying a simple cut at Mb — 5 log h — 
—20.7 + Q(z — 1.0). Table [T] summarizes the properties of 
the resulting samples, and Figure [3] shows colour-magnitude 
and colour-redshift diagrams for the galaxies in each sample. 
These samples should also be useful for future studies of 
colour evolution in DEEP2. 

Selecting galaxy samples according to Equation [1] ef- 
fectively creates a sample that is volume-limited for each 
colour. That is, the selection is colour-dependent, but the 
magnitude limit relative to M* is redshift-independent. This 
selection is different than the traditional volume-limited se- 
lection, but it nevertheless allows for comparison of the rela- 
tive numbers of galaxies of different colours at different red- 
shifts; that is, it allows us to examine the evolution of the 
blue fraction. It should be noted, however, that our values 
of fb cannot be meaningfully compared to published values 
at low redshift, which have typically been computed using 
a simple magnitude cut. In particular, our values of will 
be much higher than local values from the literature since 
the magnitude cut in Equation [T] selects many more blue 
galaxies than red. Indeed, fb values cannot even be com- 
pared among the different samples defined in Table [1] The 
absolute values of fb in each sample are mainly set by the 
selection; useful information comes only from splitting these 
samples into subsamples to look for trends with redshift or 
with galaxy properties. 

Before we move on, it is important to note that the ap- 
plication of our colour-magnitude selection criteria has im- 
portant implications for the definitions of our group and field 
samples. The VDM group catalogue has been defined using 
the entire apparent-magnitude- limited survey. It is thus pos- 
sible that a group could have been detected with two or more 
members at z = 0.8, say, while only one of these members is 
brigher than M cut . Had this group instead been at z — 1.3, 
it clearly would not have been identified, and its one bright 
galaxy would have been considered to lie in the field. If no 
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Figure 2. Rest-frame colour-magnitude diagrams for DEEP2 galaxies in redshift bins of width 8z = 0.05. Crosses (red) denote galaxies 
in groups, and diamonds (blue) show field galaxies. Groups that have fewer than two members above the magnitude limit have not been 
excluded from the group sample in this plot, although they are excluded elsewhere in the paper, as discussed in the text. The dashed 
lines show the effective DEEP2 apparent magnitude limit, as given in Equation[T] using parameters from Table[T]for Zn m = 1.3. Dotted 
lines indicate the division between blue and red galaxies as given in Equation [5] 



correction is made for this effect, the application of M cu t will 
cause a redshift nonuniformity in our samples, with galaxies 
from relatively faint groups being counted as group mem- 
bers at low z and as field members at high z. To avoid such 
biases, we redefine the group sample to be only those galax- 
ies that reside in groups with two or more members brighter 
than Mcut. The field sample comprises all other galaxies. 

We also define a group's richness N to be the number 



of galaxies with Mb ^ M cut instead of the total number of 
observed galaxies, so that richness is measured consistently 
at all z (other group properties, like a v and mean z, are still 
computed using all observed group members, for the sake of 
robustness). This group definition ensures that our sample 
of group galaxies is drawn from comparable groups at all z, 
although it means that a given group may have a different 
value of N in each sample. Also, throughout this paper we 
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Figure 3. Summary of the galaxy catalogue in each of the four 
samples defined in Table [T] The left-hand panels show colour- 
magnitude diagrams for the samples, and the right-hand panels 
show colour-redshift plots. The redshift limits of the samples are 
apparent, as are the effects of our colour-magnitude selection cri- 
teria (Equation [TJ . In the bottom three rows, contour lines are 
evenly spaced in point density. 



will restrict the group galaxy sample to those galaxies whose 
host groups have velocity dispersions a v ^ 10 km s" 1 . This 
is because, as shown in iGerke et all (|2005l ) (see figure a 
of that paper), groups with lower velocity dispersion are 
predominantly false detections. We therefore class galaxies 
in such groups with the field sample. This point is discussed 
further in § 15.11 Figure [4] shows the distribution of groups 
and their member galaxies as a function of N and a v for 
each of the four samples. 

4.2 Correcting for incompleteness 

Before we may proceed with measuring fb, an additional ef- 
fect must be accounted for. Because of the finite slit length 
available on the DEIMOS spectrograph, DEEP2 can target 
for spectroscopy only ~ 60% of the galaxies that meet its 
selection criteria (we will call the remaining potential tar- 
gets "unobserved galaxies"). Moreover, ~ 30% of galaxies 
that are targeted for spectroscopic observation fail to yield 
redshifts (we will call these "redshift failures"). Follow-up 
observations have shown that many (~ 50%) of the redshift 
failures, especially for blue galaxies, lie at redshifts beyond 
the range probed by DEEP2 (C. Steidel, private communica- 
tion). The overall failure rate is also correlated with observed 
galaxy colour and magnitude, although the trends are only 
slight. Other failures may occur because of poor observing 



Figure 4. Summary of the group catalogues in each of the four 
samples defined in Table [T] Solid lines show the differential dis- 
tribution of groups as a function of richness (number of galaxies 
above M cu t for the sample) and velocity dispersion a v . Dashed 
lines show the total numbers of galaxies within those groups that 
fall in a given bin. 



conditions, data reduction errors (e.g., poor subtraction of 
night-sky emission), or instrumental effects (e.g., bad CCD 
pixels), but such problems are only a minor contribution to 
the overall failure rate. These various sources of incomplete- 
ness may introduce spurious evolutionary trends in fb unless 
an appropriate correction is applied. 

In this work we adopt t he corrective weighting scheme 
used bv lWillmer et al.l l)200rj ) in measuring the DEEP 2 lumi- 
no sity functi o n. A m ethod of this sort was first implemented 
bv lLin et al.1 [|l999h for the CNOC2 Redshift survey. We lo- 
cate each galaxy in the three-dimensional space defined by 
observed R magnitude and R—I and B — R colours. We then 
define a cubical bin (0.25 magnitudes on a side) in this space 
around each galaxy with a successful redshift and compute 
a weight for the ith such galaxy: 



Xi = 1 + 



J2j P i( Z lo < Zj < 2hi) 



Ni 



(2) 



Here Nl is the number of galaxies in the bin around galaxy 
i with successful redshifts within the nominal DEEP2 red- 
shift range {z\ a — 0.7 ^ z ^ 1.4 = Zbi), j is an index that 
runs over all galaxies in the bin for which DEEP2 did not 
successfully measure a redshift (whether they were observed 
or not), and Pi(z\ ^ Zj ^ Zhi) denotes the probability that 
galaxy j (in the bin around galaxy i) falls in the nominal 
DEEP2 redshift range. 

To compute the probability P in equation [2] we must 
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Table 1. Summary of the data samples 





Sample I 


Sample II 


Sample III 


Sample IV 


Description 


Simple M B - M* 
limit to z=1.0 


Colour- magnitude 
complete to z=1.0 


Colour-magnitude 
complete to z=1.15 


Colour-magnitude 
complete to z=1.3 


a 
b 
c 
d 


1.0 


-20.70 


1.0 
-1.34 
-18.55 
-2.08 
-17.75 


1.15 
-1.55 
-18.77 
-2.32 
-18.16 


1.3 
-1.94 
-18.92 
-2.90 
-18.53 


# galaxies 

# groups^ 

# group galaxies 


2691 
232 
531 


9546 
863 
2588 


11767 
933 
2605 


12493 
851 
2211 


overall f[F 
field f b 
group /(, 


0.603 ± 0.011 
0.624 ± 0.012 
0.517 ± 0.025 


0.827 ± 0.005 
0.854 ± 0.005 
0.749 ± 0.010 


0.877 ± 0.004 
0.893 ± 0.004 
0.818 ± 0.008 


0.924 ± 0.003 
0.933 ± 0.003 
0.876 ± 0.008 



a The parameters 2rj m , a, 6, c, and d are defined in Equation [T] 

^ Groups must have two or more members above the sample's magnitude limit. 

c See Equation [6] for the definition of the blue fraction, /;,. 



construct a model for the redshift distributions of redshift 
failures and of unobserved galaxies. Modeling this in detail 
would require additional assumptions, but it is reasonably 
certain that the truth lies between two extreme models. As 
noted previously, many of the failures are known to be at 
z > 1.4, so we can start by making the extreme assumption 
that all failures lie in this range. In this case, called the 
"minimal" model, all redshift failures have Pi(z\ ^ Zj ^ 
«hi) — in equation [2] and unobserved galaxies have 

Pi{zi ° * Zj < Zhi) = m + N :i + Ni h + Ny (3) 

where N l zl is the number of galaxies (in the ith bin in colour- 
colour-magnitude space) with redshifts observed to be below 
zio, Nl h is the number with z > Zhi, and Nl is the number 
of observed galaxies that failed to yield a redshift. Taking 
the opposite extreme, we could assume that redshift failures 
have exactly the same redshift range as the galaxies with 
successful redshifts. In this case, called the "average" model, 
both redshift failures and unobserved galaxies have the same 
value for Pf. 

Pi(z\o ^ Zi ^ Zhi) = i ^ — • (4) 

Since blue galaxies with failed redshifts ar e known to be 
freque ntly beyond the DEEP2 redshift range. Iwillmer et al.l 
(2006) adopted a compromise model (the "optimal" model), 
in which blue galaxies (see § 14.31 for the precise definition) 
are corrected with the minimal model while red galaxies 
are weighted using the average model. We shall adopt this 
scheme in what follows; however, our basic conclusions are 
insensitive to the weighting scheme used and in fact are un- 
changed even when no weighting is used. 



4.3 Computing the blue fraction 

Having defined a galaxy sample using equation[T]and a set of 
galaxy weights \i using equation [2l we may now c ompute ft 
for our chosen sample. As in W illmer et al.l {200(f) we divide 
the galaxies into red and blue subsamples according to the 
observed bimodality in galaxy colours. This division corre- 
sponds to the dotted line shown in all panels of Figure [5] 
given by the equation 

U-B = -0.032(Ms - 5 log h + 21.62) + 1.285 - 0.25, (5) 

which was derived from the Ivan Dokkum et"aH (|2000h 
colour-magnitude relation for red galaxies in distant clus- 
ters. Eqn. [5] shifts that relation downward by 0.25 mag- 
nitudes to pa ss through the valley in the colour distribu- 
tion, following I Willmer et al.l (2006). One can also allow this 
colour division to evolve with redshift according to passive 
stellar evolution models; however, there is no evidence in the 
DEEP2 data for any evolution in the position of the valley 
(See Figure [2}, and a realistic amount of evolution has only 
minimal effects on our results (see §5.2). In the interest of 
simplicity, therefore, we do not allow the division to evolve 
in most of what follows. This division between red and blue 
galaxies produces a set of M blue galaxies out of a total 
sample of iV galaxies, each of which is assigned a weight \i- 
The corrected blue fraction is then given by 

ft = Aft /Aft, (6) 

where the number of blue galaxies Aft and the total number 
of galaxies Aft are defined as 

M N 
M b = ^2xj,Mt = ^ X i (7) 

3=1 i=l 

with the index j running over all of the blue galaxies and 
the index i running over the full galaxy sample. 
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4.4 Estimating errors on ft 

The formal error in fb is given by simple binomial statistics: 
2 N b (N t -N b ) 

^bin(fb) - ^=3 (°J 

(except in the case Afb = 0, for which |5e Propris et al.l 
(2004) argue that a(ft) = l/(2A/t))- However, this formula 
will not fully account for the scatter in our measured values 
of ft. Because a correlation exists in DEEP2 between lo - 
cal galaxy density and galaxy colour l|Cooper et alj|2006bl ). 
large-scale structure will induce an intrinsic scatter in val- 
ues of ft measured over a finite volume, in addition to the 
formal binomial error; this can be thought of, in essence, as 
a contribution to the error from cosmic variance. We discuss 
the effects of this scatter further in § 15.21 Also, errors in the 
galaxy weights %i wm a dd scatter to the measured /(, values. 
We therefore find it most convenient to estimate our errors 
empirically. 

We do this in three different ways to ensure that our 
methods are robust. First, we compute fb for each of the 
10 DEEP2 photometric pointings individually and average 
those computed values; the error o~(fb) is then taken to be 
the standard error on that mean. Second, we estimate cr(fb) 
using a jackknife sampling strategy in which each pointing 
is removed in turn from the full dataset being considered, f b 
is computed for each subsample, and the error on for the 
full sample is estimated using the usual jackknife error for- 
mula. Finally, we compute the errors on fb using a bootstrap 
resampling strategy in which the entire data subsample un- 
der consideration (e.g. all group galaxies in a given redshift 
bin) is randomly sampled with replacement, and the error 
on fb is estimated from the standard deviation of the fb dis- 
tribution for 500 such Monte-Carlo samples, using standard 
bootstrap methods. In each case, the estimated error must 
be scaled up by a factor of 1.08 to account for the covariance 
that exists between contiguous pointings due to large-scale 
structure fluctuations (i.e., cosmic variance). This factor is 
derived from Monte Carlo tests based on the covariance be- 
tween fields with the actual DEE P2 geometry, using the c os- 
mic variance calculation code of lNewman fc David l|2002l ): it 
corresponds to the conservative assumption that pointing- 
to-pointing variations in fb are dominated by cosmic vari- 
ance. 

These three different methods give comparable results, 
though the standard-error and jackknife estimates are much 
noisier than those from bootstrapping. In the interest of sta- 
bility, we will always report the bootstrap error values in 
what follows. It is, however, worth emphasizing that even 
this method may not account for all sources of fb variance 
in a given redshift bin, since the DEEP2 dataset is finite. 
Because fb seems to be lower in higher- density regions, if 
we are particularly unlucky and there is a net overdensity 
or underdensity at a given z in most of the DEEP2 fields, 
this may lead to a fluctuation in fb that is not accounted for 
in our error bars. In particular, it appears that there may 
be an underdensity at z ~ 0.9 in most of the DEEP2 fields; 
this may lead to unusually high values of fb at that redshift. 
To ameliorate this problem, in addition to computing fb in 
independent bins, we will also sometimes compute fb in a 
sliding box, which will smooth out any large-scale structure 



fluctuations at the expense of introducing bin-to-bin corre- 
lations in the resulting measurements. 



5 RESULTS 

Table [1] shows the properties of the four samples defined in 
§ 14.11 including blue fractions. A primary result of this work 
is apparent in the last two rows of the Table: in all four 
samples considered, the blue fraction is significantly lower 
in groups than it is in the field. Thus, the well-known quali- 
tative distinction between local field and group (or cluster) 
galaxy populations was in place by z ~ 1. It is, however, 
worth emphasizing that the quantitative increase in val- 
ues from sample I to sample IV is not evidence for evolution 
because the overall value of in each sample is strongly 
affected by the sample-selection procedures in § 14.11 (i.e., 
sample IV has a much more sharply tilted selection cut than 
sample II, so it will include proportionally fewer red galax- 
ies by construction). In what follows, we will divide each 
sample into redshift bins, allowing for a meaningful inves- 
tigation of evolutionary trends in the blue fraction. First, 
though, it is important to understand the trends that exist 
at fixed redshift. 

5.1 Blue- fraction trends at fixed redshift 

As discussed in the introduction, the blue fraction in local 
clusters exhibits a large system-to-system scatter that arises, 
at least in part, because the measured value of ft depends on 
exact ly which galaxies and clusters are used to co mpute it 
(e.g., |Pe Propris et alj|2004l : iPoggianti et~al"1l2006T l. Before 
considering fb evolution in DEEP2, then, it is crucial to 
understand trends in the blue fraction at fixed redshift; oth- 
erwise they may introduce selection effe cts that enhance o r 
compete with evolutionary effects (e.g.. Ismail et af]ll998l ). 
As discussed in § !5.2l below. there is little to no evolution in 
fb over the redshift range 0.75 < z < 1.0 in DEEP2, so we 
will limit ourselves to this range (i.e. to samples I and II) 
when studying trends at fixed redshift. 

First, it is interesting to consider the effect of galaxy 
magnitude on the blue fraction. In doing this, it is important 
to use the same absolute magnitude cuts at all colours (i.e., 
vertical lines in Figure [2}. If one instead used tilted cuts like 
the ones that define samples II-TV (Equation[T} , the extreme 
slopes of these cuts in colour-magnitude space would mean 
that the brightest magnitude bins would exclude all red 
galaxies while keeping some blue ones, leading to obviously 
spurious results. Thus we are limited here to subsamples of 
sample I, the smallest sample in this study. It is also impor- 
tant to allow the magnitude bins to evolve with redshift in 
the same manner as M* , to ensure that similar galaxies are 
being compared across the redshift range. Thus, we divide 
sample I into bins in the quantity Mb — 5 log h + Q(z — 1), 
with Q = -1.37. 

Figure [5] shows ft in bins of absolute magnitude for 
galaxies in groups and in the field. Groups here are defined to 
be those systems that have two or more members in sample 
I. A trend is evident for both populations, with brighter sub- 
samples having a lower blue fraction — except for the bright- 
est galaxies, where the trend may reverse. T h is last result 
is consistent with the results of lCooper et al.l (|200 6b). who 
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Figure 5. Dependence of the blue fraction fi on absolute mag- 
nitude in sample I (0.75 ^ z ^ 1-0), for galaxies in groups (dia- 
monds with horizontal bars) and in the field (triangles). The data 
points show in bins of evolution corrected absolute magnitude, 
Mb — 5 log h + Q (z — 1.0), where we have adopted the lFaber et al.l 
(2006) value of Q = —1.37. Horizontal bars (suppressed for the 
field sample) show the bin size used to compute each data point, 
and vertical bars show the la error on fb, computed from boot- 
strap resampling of the galaxies in each bin. Data points have 
been offset slightly in the horizontal direction for clarity. Groups 
contain fewer blue galaxies than the field at all magnitudes except 
possibly in the brightest bin. In addition, a trend is apparent for 
both group and field galaxies, with brighter galaxies being redder 
on average, except at the bright end, where the trend may reverse. 
Linear fits to both trends, excluding the brightest bin, show that 
their slopes are consistent with the same value. 



find a population of bright, blue DEEP2 galaxies in high- 
density environments; it can possibly be ascribed to AGN 
and starburst activity. It is also interesting to compare the 
trends for the fainter group and field populations. Fitting 
a straight line for each subsample, excluding the brightest 
bin, yields a slope of 0.102 ± 0.033 for the field population 
and 0.106 ± 0.044 for the group galaxies. That is, there is 
evidence for a similar trend with limiting magnitude for /;,, 
both in groups and in the field. We will discuss the implica- 
tions of this result in § [7j 

We will use Sample II to explore other possible system- 
atic trends in at fixed z, since it is significantly larger than 
sample I by virtue of the tilted colour-magnitude selection 
criterion used to define it. Figure [6] shows fb for the group 
galaxy population only, binnned by various group or galaxy 
properties. For comparison, the blue fraction for field galax- 
ies is shown as a dotted line. The upper left panel shows the 
dependence of fb on the velocity dispersion a v of the galax- 
ies' host groups — that is, each bin contains all of the galaxies 
in groups with a v in the range of that bin. In this panel, we 
have included in the group sample those galaxies in groups 
with a v < 100 km s _1 , unlike all other plots in this paper. 
As shown, in terms of ft these galaxies bear a much stronger 



resemblance to the field than to the rest of the group sam- 
ple. This is no surprise, since these groups are known to be 
predominantly false detections — i.e., chance associations of 
field galaxieaj (|Gerke et alj 1200*51 '). Figure [6] then stands as 
further justification for our decision to class these galaxies 
with the field. The other data points in this panel show an 
interesting trend with a v : fb declinines with a v at low dis- 
persion and then rises again at high dispersion. Caution is 
warranted, however, because (1) high dispersion groups are 
more likely to be contaminated by interloper field galaxies, 
since the VDM group finder uses a larger search volume for 
such groups; and (2) the large uncertainty in the measured 
values of a v for DEEP2 groups makes it difficult to assess the 
reality of observed trends of galaxy properties with group ve- 
locity dispersion. In particular, point (2) might be a problem 
because the sample has many groups with two members, for 
which measurements of a v are maximally uncertain. Such 
groups might dominate over the population of true high- 
dispersion systems, which are rare. If there were, say, an 
underlying monotonic decline in fb with increasing a v , then 
the up-scattered small groups could induce an apparent up- 
turn at the high-a„ end. For these reasons, we refrain from 
drawing conclusions regarding the relation between fb and 
a v . 

There is, however, an apparent trend with group rich- 
ness, iV, at high values of N, visible in the top right panel of 
the figure. Richness in this context is defined to be the num- 
ber of galaxies in a given group above the colour-magnitude 
limit M cu t that defines sample II (see Equation [1] and Ta- 
ble [TJ. The blue fraction declines at the highest richness 
values. This result illustrates the primary reason that we 
included in our sample only those groups with two or more 
members above M cu t: had we not done so, we would have 
been sampling richer groups at higher redshifts, potentially 
causing a spurious decrease in fb with increasing redshift. 

The remaining four panels investigate the dependence 
of fb on group members' distance from their host groups' 
centres (group-centric radius) and on their peculiar velocity 
relative to their host groups. To probe dependence on group- 
centric radius, we compute fb for group galaxies within an- 
nuli on the sky around the mean right ascension and dec- 
lination of their host group. The annuli are defined in two 
ways — in units of comoving Mpc and also in units of r2oo, 
the radius at which the group is 200 times denser than the 
background. This radius can be estimated from the group's 
radial velocity dispersion ay as 



T"200 = 



\/3cr r 
WH{z) 



(9) 



where H(z) is the r edshift-dependent Hubble parameter 
l|Carlberg et alJll997h . We compute each group member's 
peculiar velocity with respect to the mean redshift of its 
host group as v pcc = c(z — z)/(l + z), and we consider the 
peculiar velocity both in units of km s _1 and normalized to 
the host group's velocity dispersion a v . The middle two pan- 
els of Figure [6] show the dependence of fb in groups on the 
group-centric distance, and the bottom two panels show the 



2 it is true, though, that such extremely poor groups in 
the local Universe are in fact frequently dominated by spi- 



rals jZabludoff fc Mulchaevlll99d) . 
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Figure 6. Dependence of the blue fraction /(, of group galaxies in sample II (0.75 ^ z ^ 1.0) on various group parameters: from upper 
left to lower right, these are the velocity dispersion of the host group, cr v \ the richness (number of galaxies above the magnitude limit) of 
the host group; distance from the centre of the host group, in units of T200 (see text) and in megaparsecs; and peculiar velocity relative 
to the host group's redshift, in units of the host group's velocity dispersion, and in km s _1 . In each panel, the data points show /j, in 
bins whose sizes vary to ensure robust statistics; the bin sizes are shown by the horizontal bars. As in Figure [5] vertical bars show the 
1<t bootstrap errors. The value of for field galaxies is shown as the dotted line. See § I5.1l for details of the selection in each panel and 
§ 14.31 for details of the computation of ft,. 



dependence on peculiar velocity. None of these plots shows 
a significant trend, and linear fits to the data points in each 
panel are consistent with zero slope. 

The lack of a trend with group-centric radius may 
be surpri sing in light of studies of groups and cluste rs at 



low (e.Q.. IWhitmore et al 
high (e.g., iPostman et al 



1995j lDe Propris et al.|[200j ) and 



20051 ) redshift, showing a strong 



relation between cluster-centric radius and galaxy type. 
However, there is a large uncertainty in the determinations 
of DEEP2 group centres and mean redshifts since these in- 
volve taking means of a small number of objects, so it is 
likely that any trends would be significantly diluted by scat- 
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ter in these values[f] Nevertheless, there is a hint in the lower 
left panel that galaxies at low values of v pec /(T v and r/r2oo 
tend to be redder than the general group population. In ad- 
dition, in the lower right panel it appears that group galaxies 
with moderate absolute peculiar velocities are redder than 
usual, while those with the highest velocities are blue. This 
echoes the trend with a v at upper left (as expected, since 
only high-crt, groups contain galaxies with high values of 
Wpec); these two panels taken together suggest that a nega- 
tive correlation between fb and a v for DEEP2 groups may 
be masked by the presence of interloper field galaxies that 
preferentially lie at high values of v pec . It would be preferable 
if we could compute Fig. [6] using only high-richness groups, 
whose dispersions and centers can be determined more ro- 
bustly; however, restricting even to groups with N > 3 re- 
duces the sample size so significantly that no firm conclu- 
sions can be drawn in light of the resulting error bars. Thus, 
in the absence of any significant measured trends, we see 
no justification for further restricting our sample by pecu- 
liar velocity or group-centric radius — for instance by taking 
only galaxies within r2oo as has been done in various othe r 
studies (e.g., |Pe Propris et al.ll2004l ; iPoggianti et al.l l2006). 
in which it was possible to measure r2oo more robustly than 
it is here. In what follows, the sample of group galaxies will 
always include all galaxies in groups with a v > 100 km s _1 
and with two or more members above M cut , regardless of 
radius or peculiar velocity. 

5.2 The evolution of the blue fraction 

Having ensured that the group sample will not be affected 
by trends at fixed redshift that will complicate comparisons 
at different redshifts, we now proceed to probe the evolution 
of the blue fraction with z. We divide each of the samples 
in Table [T] into several independent redshift bins containing 
equal numbers of galaxies, and we compute ft in each bin 
for the group and field galaxy subsamples. As discussed in 
§ 14.41 there may be scatter from bin to bin that exceeds 
the error estimates, especially at z ~ 0.9, so to smooth out 
this effect we also compute fa in a sliding box whose width 
is twice the average width of the independent bins in each 
sample. 

The resulting fb values for each sample are shown in 
Figure [7] as a function of redshift and of cosmic time; this 
Figure summarizes the principal results of this paper. As 
noted previously, there is clearly a significant difference in 
the blue fractions of group and field populations at z ~ 
1. In addition, in all four samples there is no evidence for 
evolution in fb over the range 0.75 < z < 1.0: group and field 
populations are consistent with constant ft over this range. 
However, the bottom two panels of the Figure show that fb 
evolves dramatically at z > 1.0, and this evolution is much 
stronger in groups than it is in the field. In fact, interestingly, 
the two populations appear to converge in terms of fb at 
z ~ 1.3. It is also notable that the evolution in samples 

3 Indeed, the particularly attentive reader will have noticed that 
the highest measured values of r/r200 here are ~ 10, which is ex- 
tremely large for any realistic system; such values are attributable 
to groups with spuriously low measured values of u v and hence 
f200- 



III and IV appears to flatten below z ~ 1, in agreement 
with samples I and II. It is important to emphasize that the 
bins used in the lower right panel of Figure [7] are broader 
than those shown in Figure [2] so that each bin contains a 
substantial number of red galaxies; that is, the convergence 
of the group and field populations at z ~ 1.3 does not result 
from small-number statistics. 

To further assess the robustness of these results, we have 
investigated the effects of changing our assumptions about 
the DEEP2 redshift incompleteness modeling and the lumi- 
nosity evolution of the sample. In each panel of Figure [7J 
th e galaxy weigh t ing sc heme used is the "optimal" scheme 
of IWillmer et al] (|2006T ). and the mag nitude evolu ti on pa - 
rameter is set to the value found in iFab cr et all (2006), 
Q = —1.37. Changing the w eighting scheme to eith er the av- 
erage or minimal models of lWillmer et all (I2006T ) — or even 
using no weights whatsoever — effects no qualitative change 
in the results shown in Figure [7] Similarly, changing the 
value of Q by ±0.31 (the ler uncertainty on this parame- 
ter from iFaber et al.l 120061 ) has no qualitative effect on the 
results. 

As a further test, we also consider the possibility of 
evolution in the colour of the "green valley" separating red 
and blue galaxies in colour-magnitude space. If the colour of 
the valley changes with time, the fb evolution seen in Fig- 
ure [7] might simply be the result of passive evolution mov- 
ing galaxies redwar d across our (fixed) dividing line (e.g., 
lAndreon et al.ll2006t ). However, the existence of a bimodality 
in the galaxy colour distribution suggests that the division 
between red and blue galaxies should be drawn through the 
valley, as we have drawn it, rather than being tied to the 
colours of red galaxies alone. Simple inspection of Figure [2] 
shows no clear evidence for evolution in the locus of this 
valley. Indeed, if one adopts a model in which red galaxies 
are continuously being formed from blue ones throughout 
the redsh ift range of interest [e.g., the "hybrid" model pro- 
posed bv lFaber etail (|2006T )1. then one would expect little 
evolution in the valley's position. As old red galaxies evolve 
passively toward redder colours, new red galaxies are formed 
to take their place, resulting in a red sequence that broad- 
ens with time, while the valley between red and blue galaxies 
remains nearly fixed. 

Nevertheless, it is wort h exploring the effect of allowing 
the dividing line to evolve. iBlantonl (|2006h showed that the 
locus of the green valley has evolved redward by at most 
~ 0.1 magnitudes in u — g from z = 1 to the present day. 
We test the effect of this evolution by applying a colour 
separation which is given by Equation [S] at z = 1, with a 
linear evolution of 0.1 magnitudes redward per unit decrease 
in redshift (that is, the blue-red dividing line gets redder 
with time, as might be expected from passive evolution). The 
evolution of fb with this evolving cut in sample IV is shown 
in Fig [8] As shown, the evolving colour separation produces 
a (marginally significant) drop in the group fb from z — 0.75 
to z = 1.0, probably because the green-valley evolution we 
assume here is too strong. This is followed by a rise at higher 
redshifts and a convergence of the group and field values at 
z ~ 1.3. Hence, the main results of Figure [7] still hold: there 
is substantial evolution in fb in groups at z > 1.0, with a 
convergence of the group and field fb values at z ~ 1.3. 

Before this result can be fully accepted, however, there 
is another systematic effect that should be considered. It 
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Figure 7. Evolution of the blue fraction with redshift and cosmic time (time since the big bang). For all four samples defined in Table[T] 
the evolution of is shown from z = 0.75 to the limiting redshift of each sample. In each panel, data points show values in six 
independent bins (only four for sample I, which is much smaller than the others) containing equal numbers of galaxies (Note that these 
bins are typically broader than the ones used in Figure [3] each bin contains a significant number of red galaxies) . Blue triangles show /(, 
for field galaxies, and red diamonds denote for group galaxies; these points have been slightly offset in the redshift direction for clarity. 
Lines show lc errors on computed in a sliding bin of twice the mean width of the independent bins. Note that the four panels' axes 
have different horizontal and vertical scales; the absence of the upturn at high redshift in the top two panels is owing to the fact that 
these panels do not extend beyond 2 = 1. Over the redshift ranges in common, all samples are fairly consistent. In particular, no sample 
shows significant evolution in over the range 0.75 < z < 1, but for 2^1 (samples III and especially IV), /(, evolves substantially, with 
much stronger evolution in groups than in the field. The group and field samples become nearly indistinguishable at the highest redshifts 
considered. See § 14. 3 1 for details of the computation of the blue fraction. 



has been observed (e.g, iBell et alj l2004al ; IWeinerl 120051 ) 
that there are two distinct classes of red galaxies at high 
z: nonstar-forming, "red-and-dead" early-types and dust- 
reddened, star-forming late-types. The latter class will pre- 
sumably exhibit spectral emission lines; hence it should be 
relatively easy to obtain redshifts for these galaxies, even if 
their spectra have low signal-to-noise ratios. Redshifts may 
be harder to obtain for the nonstar-forming class, since these 
exhibit weaker spectral features in generafl. It is thus pos- 



sible that DEEP2 fails to obtain redshifts preferentially for 
high-redshift (i.e. faint) early-types. This would mimic evo- 
lution in both in groups and in the field, but if early- 
type galaxies occur preferentially in groups, this selection 
effect could potentially also lead to the convergence seen in 
the group and field /j, values at z ~ 1.3. In principle, the 
weighting scheme described in § 14.21 should approximately 
correct for this effect; however, if the situation is so dire 
that all galaxies at a given redshift and in a particular re- 



4 However, it has recently been shown llYan et alj|200rj) that ~ 
40% of nonstar-forming galaxies in the SDSS exhibit significant 



[O II] emission from AGN activity; redshifts should be readily 
obtainable for these galaxies. 
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Figure 8. Evolution of /{, in sample IV, when the dividing line 
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0.1 magnitudes per unit decrease in z I IBlantonl 120061) . There is 
a (probably spurious) drop in /j, out to z = 1.0 (see text), but 
at higher redshifts the evolution is the same as observed with a 
non-evolving separation in Figure [7] 



gion of colour-colour-magnitude space are redshift failures, 
then the weighting will be of no help. 

This potential bias must therefore be taken seriously. 
Fortunately, it is possible to get some sense of its impor- 
tance by looking at the HST/ ACS imagi ng that the AEGI S 
collaboration has obtained in the EGS l|Davis et alj|2006h . 
Supposing that the true mix of spheroidal and dusty-spiral 
red galaxies in the field is constant with redshift, if early- 
types are being preferentially missed at high redshift, this 
would result in a declining ratio of dead to dusty among red 
field galaxies in DEEP2. We have examined all of the red 
field galaxies with HST imaging and confirmed DEEP2 red- 
shifts in the range 0.75 < z < 1.3 in the EGS (a total of 63 
galaxies). At all redshifts, roughly 50% of these show no vis- 
ible evidence for star formation or spiral structure. Hence, 
it appears that DEEP2 is not preferentially failing to ob- 
tain redshifts for a particular class of red galaxies at high 
z, although it is worth emphasizing that this conclusion has 
been drawn from a limited sample of galaxie^f]. 

The strong evolution of the DEEP2 group blue fraction 
at z > 1.0 therefore appears to be robust. It can be explained 
either by a blue galaxy population that decreases with time, 
a red galaxy population that in c reases with time, or both. 
As shown in iBell et al.l (|2004bl) . IWillmer et all (|2006l ) and 
iFaber et al.1 ( 20061 ). there is significant growth in the typ- 
ical comoving number density <j>* of red-sequence galaxies 
over the redshift range probed by DEEP2, while (j)* for 
blue galaxies has remained roughly constant over this range. 



5 It may be possible to address this issue more directly in the 
near future, by using Spitzer /IRAC data to obtain photometric 
redshift information for DEEP2 redshift failures. 
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Figure 9. Incompleteness-corrected galaxy counts used to com- 
pute for sample IV in Figure [7] Solid lines show the total 
number of galaxies, Aft (see Equation (7)1 , in groups (top) and the 
field (bottom), for the six redshift bins shown in the bottom right 
panel of Figure [7] Dotted lines show the number of blue galax- 
ies, Afb, in these bins, and dashed lines show the number of red 
galaxies, Af r = Aft —Aft- As shown, while the overall number of 
group galaxies is growing with time, the number of red galaxies 
in groups increases much more rapidly, increasing by almost an 
order of magnitude over the redshift range shown. 



Hence, it seems reasonable to ascribe the evolution of ft 
to the preferential build-up of the red galaxy population 
in groups. This conclusion is made manifestly clear in Fig- 
ure [9] Rather than showing fractions, this figure shows the 
(incompleteness-weighted) numbers of galaxies from sample 
IV in each of the redshift bins shown in the lower right panel 
of Figure[7] both in groups and in the field. More specifically, 
the figure shows the total number of galaxies Aft, the number 
of blue galaxies Afb (see Equation 0, and the number of red 
galaxies Af r = Aft — Afb, all as a function of redshift, both in 
groups and in the field. It is important to reiterate that the 
redshift bins were chosen to contain equal unweighted total 
numbers of galaxies, so the overall weighted numbers in each 
bin are not constant. In any case, we find that the results of 
this paper are insensitive to the weighting scheme we choose 
and remain unchanged even if no weights are applied at all. 

The figure shows that Aft for groups increases with time, 
as expected in the hierarchical ACDM structure-formation 
paradigm. But Af r in groups grows at a much faster rate, 
increasing by almost an order of magnitude from z ~ 1.3 
down to z ~ 0.75. At the same time Aft in the field stays 
nearly constant, and Af r in the field increases only modestly 
(by a factor of < 2) over this range. It seems clear, then, 
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that the build-up of red sequence galaxies over the DEEP2 
redshift range has taken place preferentially within groups 
and clusters. 



6 TESTS WITH MOCK CATALOGUES 

The above results so far appear to be robust, but to es- 
tablish them definitively it is vital that we also test for 
biases introduced by the group-finding procedure. For ex- 
ample, it is well known (and has been confirmed here) that 
groups and clusters are populated by redder galaxies than 
the field, but, as shown in Figure [2] the DEEP2 apparent 
magnitude limit selects strongly against red galaxies at high 
redshift. If red galaxies dominate groups and clusters, this 
effect might cause their richnesses to appear strongly re- 
duced at high redshift — so much so that they might eventu- 
ally become undetectable with the VDM group-finder. The 
high-redshift VDM group sample would then be dominated 
by false detections. If this were the case, then the group and 
field populations identified by the VDM would necessarily 
be indistinguishable, since they would essentially be random 
subsamples of the same galaxy population, and so the con- 
vergence of the group and field populations seen in Figure [7] 
would simply be the result of group-finding errors and not 
of genuine galaxy evolution in groups. 

We can test for such an effect by running the VDM 
group-finder on the DEEP2 mock catalogues described in 
§ [3] As discussed in that section, the mock galaxies are as- 
signed rest-frame U — B colours according to the observed 
DEEP2 colour-magnitude-environment in the redshift range 
0.8 < z < 1.0. Figure [101 shows the evolution of fb in the 
mocks for galaxies in real groups (i.e., galaxies that actu- 
ally reside in dark matter haloes that contain other galaxies, 
lower curve), for galaxies in found groups (i.e., galaxies in 
groups identified by the VDM algorithm, diamond points), 
and for the field samples that complement each group sam- 
ple (upper curve and triangle points). In computing fb, the 
same limit in colour-magnitude space is applied as the one 
used for sample IV (given by Equation [JJ and Table [1] and 
shown in Figure(2]|, and, as in the data, galaxies are included 
in the group sample only if their host group has two galaxies 
brigher than the limit. As Figure [10] illustrates, the VDM 
group-finder reduces the separation in between the group 
and field samples (as expected, since the found groups are 
contaminated by interloper field galaxies). However, it does 
not induce a spurious evolutionary trend in the blue frac- 
tion of the found group sample: the slopes of fb versus z are 
consistent for the real and found samples. Thus it does not 
appear that the VDM group finder is introducing spurious 
trends in the measurement of blue fraction evolution. 

Finally, it is interesting to note in Figure [10] that the 
values of fb for the true group and field samples in the 
mocks show clear evolution. As discussed in Section [3] we 
have not introduced any redshift dependence in the colour- 
environment relation used to assign colours to mock galax- 
ies. By design, the blue fraction will remain constant with 
redshift for comparable local density values. Therefore, the 
evolution of fb seen in the mocks can only be due to evolu- 
tion in the distribution of local densities within the mocks. 
That is, the growth of large-scale structure in the mocks 
causes the number of high-density regions to grow with time, 
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Figure 10. Evolution of the blue fraction in the mock catalogues 
described in Section [3] Solid lines show the true evolution of fb 
in the mocks, in groups (lower curve) and in the field (upper 
curve), computed using real groups, i.e., sets of galaxies that 
actually occupy common dark matter haloes. Data points show 
the reconstructed evolution of in the mocks, computed using 
found groups, i.e., groups of galaxies identified in redshift space by 
the VDM group-finding algorithm. As expected, the distinction 
between found groups and the field is smaller than the distinction 
for real groups because of interlopers; however, it is clear that the 
group-finder does not induce a reconstructed evolution that is 
stronger than the actual evolution in the mocks or as strong as 
what is observed in the data. 



so there are fewer red galaxies at high redshifts in the mocks 
simply because there are fewer galaxies in overdense environ- 
ments at high redshifts. This effect provides a partial expla- 
nation for the evolution in observed in DEEP2. However, 
because the mock group and field samples remain distin- 
guishable in terms of fb at all redshifts while the observed 
samples do not, there must also be some decrease in the blue 
fraction of DEEP2 group galaxies with time at fixed over- 
density. (This can also be seen directly in Figure[9] where the 
numbers of red galaxies in groups grow more rapidly than 
the total group population.) We will discuss this further in 
Section [7J 



7 DISCUSSION 

The main results of this paper are (1) a significantly lower 
value of fb in group versus field populations at z ~ 1; (2) 
a strong increase of fb in groups with redshift at z > 1.0; 
(3) the lack of any obvious fb evolution, for either groups 
or the field, at lower redshifts; (3) a convergence of for 
the group and field populations at z ~ 1.3; (5) a decline in 
fb in richer groups; and (6) a negative correlation between 
fb and galaxy luminosity, both in groups and in the field. 
In light of Figures [7J and [9] it appears especially clear that 
the build-up of the DEEP2 red galaxy population has taken 
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place most dramatically in groups and clusters rather than 
in the field (here and throughout, field refers to all galaxies 
not in groups). 

In interpreting this conclusion, three interrelated ques- 
tions arise. 

(i) What physical mechanisms are acting in groups and 
clusters to quench star formation and produce the evolution 
seen in § 15. 2f 

(ii) Is the evolution of fb in DEEP2 groups and clusters 
simply an extension of the similar evolution seen at low 
redshift — i.e., is this just the BO effect at high z, or are 
other mechanisms in effect at this epoch? 

(iii) If groups and clusters are the site of red galaxy for- 
mation, what is the nature of the red galaxies present in the 
field, especially those at high redshift? 

By comparing these results to other studies at high and 
low redshift — and to current theoretical models of galaxy 
evolution — we can attempt to fit our observations into a co- 
herent picture of star formation and quenching in groups 
and clusters. 

7.1 Comparison to the evolving colour-density 
relation 

As shown in § 15.21 the blue fraction in DEEP2 groups evolves 
strongly over the range 1.0 < z < 1.3, and the group and 
field populations become indistinguishable in terms of fb at 
z ~ 1.3. These results are qualitatively consistent with re- 
cent work studying the evolution of the high-redshift galaxy 
population as a function of local galaxy density, both in 
DEE P2 and in other surv eys. In a companion paper to this 
one. ICooper et all (|2006al ) show that a correlation exists be- 
tween galaxy red fraction (f r — 1 — fb) and local overdensity 
out to z > 1 but that this correlation vanishes at z ~ 1.3, 
in c lose agreement with our findings here. 

INuiiten et all (|2005h used galaxies in the photometric 
CFHT Legacy Survey to show that the fractions of red 
and morphologically early-type galaxies grow with time and 
evolve more strongly in hi g h-den sity regions than elsewhere. 
In addition, iTanaka et alj (120051 ) compared photometric ob- 
servations of two high-redshift clusters to the SDSS galaxy 
population; they find that the color-density relation has 
steepened with time and that this steepening is strongest in 
very o verdense regions. However, as shown in lCooper et al.l 
(2005), photometrically determined redshifts are not suffi- 
cient for accurate d eterminations of galaxy density. Also, 
INuiiten et al.l ((2005) consider galaxies above a set of fixed 
absolute magnitude limits rather than allowing their limits 
to evolve as M*. Since galaxies that are faint relative to 
M* tend to be bluer than average (see Figure [SJ), and be- 
cause M* is brighter back in time, a fixed cut in M includes 
more faint (blue) galaxies at high redshift, which causes fb 
to increase (spuriously) with z. This selection effect acts to 
enhance t he observed evolution . 

Also. lCucciati et al.l |2006) recently considered the frac- 
tions of galaxies in various colour ranges as a function of red- 
shift and local density in the VVDS redshift survey. They 
find that the density dependence of the blue and red frac- 
tions vanish es at z ~ 1.0, a som ewhat lower redshift than 
we find (see ICooper et al. 2006a for a detailed discussion). 
They also find that the blue fraction falls with time at all 



densities, in contrast to our finding that fb is roughly con- 
stant in the field at all redshifts. We see are two possible 
reasons for this (minor) disagreement: first, these authors 
again consider only a fixed absolute magnitude limit rather 
than taking M* evolution into account, and second, as they 
also caution, their faint, red galaxy sample is incomplete at 
the highest redshifts they consider, so any evolution of fb in 
this range will be spurious. In this paper, we have been care- 
ful to design galaxy samples that avoid such subtle selection 
effects. We are thus able to confirm that the results from 
the three papers above ( Nuiiten et al.l 120051 ; ITanaka et al.l 
120051 : ICucciati et al.ll200(: ) are qualitatively correct despite 
possibly suffering from such effects. 

We have also considered groups versus the field, rather 
than local density. This will be useful in that it allows our 
results to be interpreted in terms of the dark matter haloes 
hosting DEEP2 galaxies and compared more readily to the- 
ory (see Section \7.'A\ . Local density remains a powerful tool 
for studying galaxy populations, however, since it probes a 
wider range of environments than groups and clusters, al- 
lowing one to explore whether galaxy evolution might be 
driven by processes that occur on scales larger (or smaller) 
than the scale of groups and clusters, or whether underdense 
regions are as important to galaxy evolution as overdense re- 
gions. We therefore urge the rea der to compare our results 
to those of ICooper et al.l l|2006al ). In particular, that paper 
shows that the red fraction of galaxies in overdense regions 
declines strongly with increasing redshift, while remaining 
constant in underdense regions, and that the red fractions 
of these two environments converge at z ~ 1.3. This result is 
strikingly consistent with the lower right panel of Figure [7] 
and stands as confirmation that the results presented in this 
paper are correct. 

Bo th the overdensity considerations of ICooper et al.l 
l|2006al ) and the halo-mass considerations of this paper point 
to a partial answer to the question of which mechanisms 
are driving the evolution we observe. As discussed in § 12.21 
DEEP2 samples intermediate-mass objects (in the range 
5 x 10 12 < M/ Mq < 10 14 ), r a ther th an rich (M ~ 10 15 M G ) 
clusters. Also, ICooper et al.l (|2006al ) find that the relation 
between galaxy colour and overdensity is a smooth function 
of overdensity and does not set in only at the highest densi- 
ties. We may thus confidently conclude that the quenching 
of star formation in groups and clusters cannot be ascribed 
solely to processes that are significant only in rich clusters, 
such as ram-pressure stripping and harassment (similar con- 
clusi ons have been draw n in lower-redshift studies as well, 
e.o.. lBalogh et al.ll2002t ). 

7.2 Comparison to the Butcher-Oemler effect 

It is nevertheless tempting to identify the fb evolution seen 
in Figure [7] with a simple extension of the BO effecl|f| to 
higher redshifts. It is difficult, however, to square this in- 
terpretation with the weak or nonexistent evolution seen 



6 Recall that, in this paper, we are using the phrase "Butcher- 
Oemler effect" to refer specifically to evolution of the blue fraction 
in groups and clusters owing to the formation of large-scale struc- 
ture and the attendant lower age of typical massive haloes back 
in time (see § 0. 
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over the range 0.75 ^ z < 1.0. To be sure, this redshift 
range might not sample a long enough period of cosmic 
time to discern evolution due to the BO effect: "clas sical" 
detections of the BO effect (e.g., Butcher fc Oemlcr 
Rakos fc Schombertlll995l ; iMargoniner fc de Carvalhc 
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1984 



2000 



Kodama fc Bowerl 120011 ) have typically considered clusters 
over a redshift range < z < 0.4, which corresponds to 
3h _1 Gyr of cosmic time, whereas the range 0.75 ^ z ^ 1 
covers only 0.8/i~ Gyr. Nevertheless, at minimum, if the 
evolution observed here is simply the BO effect extended 
to high-redshift groups, the evolution becomes substantially 
more rapid at z > 1.0. 

Moreover, it is not even clear that groups of the sort 
being considered in this study would be expected to ex- 
hibit a strong BO effect at lower redshifts. As shown in 
Figure [4j the vast majority of DEEP2 groups have veloc- 
ity dispersions a v < 600 km s . Because most groups also 
have very few galaxies with spectroscopy (typically fewer 
than ten), there is also a very large scatter (typically a few 
hundred km s _1 ) in the measured dispersions. Coupled with 
the fact that the number of groups falls very steeply with 
increasing a v , this scatter implies that the measured dis- 
persions will tend to be biased high by an Eddington-type 
bias, so that DEEP2 groups are typically not as massive 
as their velocity dispersions seem to impljQ. Also, as dis- 
cussed above, our mock catalogues imply that DEEP2 does 
not sample high-mass clusters. There have been very few 
studies of galaxy evolution in the interm ediate-mass sys- 
tems we do sample, but a recent study by IPoggianti et al.l 
(2006) has shown that the fraction of [O II] emitting galaxies 
does not evolve significantly out to z ~ 0.8 for groups with 
a v < 500 kms _1 . Also, in their study o f the m orphology- 
density relation at z ~ 1, ISmith et all (|2005h show that 
galaxies in intermediate-density environments show little 
m orphological evolution down to the z ~ 0.5 epoch studied 
bv lDressler et all (|l997l ). If we naively identify intermediate- 
density environments with intermediate-mass groups like the 
ones in DEEP2, we might also not expect to find significant 
evolution in ft over the range < z < 1. Indeed, prelim- 
inary indications (Gerke et al. in prep.) are that ft is lit- 
tle different in DEEP 2 groups out to z ~ 1 than it is in 
the groups detected bv lEke et ail i|2004h in the 2dFGRS, at 
lea st for the galaxy pop ulation considered in sample I (but 
see iMartmez et~al]|2006l ) . 

In light of the theoretical basis for the BO effect, it 
is not surprising that the effect would appear weak in th e 
DEEP2 sample. As first discussed by KauH'maim (1995), 
the BO effect is a natural outcome of hierarchical structure 
formation so long as cluster environments are efficient at 
quenching star formation in galaxies. Since clusters form 
at late times, galaxies in clusters observed at high red- 
shift will have spent less time on average in the cluster 
environment than their counterparts observed at low red- 
shift. Hence, they will have had less time to feel the ef- 
fects of the physical mechanisms that are acting to quench 
star formation. This picture implies that the BO effect will 
be weaker in less massive systems, as is observed in, e.g. 



7 Note that this point does not affect the estimate that DEEP2 
groups lie in the range 5 X 10 12 < M/M Q < 10 14 ; this was derived 
from the known masses of group haloes in the mock catalogues. 



IPoggianti et all l|2006l ). since smaller systems will typically 
form earlier and then merge into larger systems, resulting a 
similarl y young galaxy populatio n at all redshifts in these 
objects. iKodama fc Bower] (|200ll ) considered a set of sim- 
ple models for the star-formation histori es of cluster galax- 
ies, and their results support the basic iKauffmannl (|l995T ) 
scheme: the observed BO effect can be explained by the de- 
creasing accretion rate onto clusters at late times, coupled 
with a declining star-formation rate in field galaxies. Nei- 
ther of these studies considered the effect of dark energy on 
this picture, but it is worth noting that, in the concordance 
ACDM cosmology, the universe becomes vacuum dominated 
around z ~ 0.4, sharply attenuating matter infall onto clus- 
ters at later times. Infall rates should be much higher at the 
redshifts considered in this paper; hence one might expect a 
weaker BO effect in DEEP 2, since the quenching blue galaxy 
population in groups and clusters is being steadily replaced 
by new infalling galaxies. 

It therefore seems reasonable to conclude that the evo- 
lution seen at z > 1 in DEEP2 is different in nature from the 
low-redshift BO effect. This is not to say that the physical 
mechanisms quenching star formation in individual galax- 
ies are necessarily different at the two epochs; rather, the 
evolution may simply be in the efficiency with which these 
mechanisms alter the bulk properties of the cluster galaxy 
population. In particular, the rapid evolution at z > 1.0 and 
the convergence of the group and field ft values at z ~ 1.3 
suggest a picture in which, at the epochs considered in this 
paper, groups have only recently become efficient environ- 
ments for quenching, while the group galaxy population has 
a similar star-formation rate to the field at earlier times. 
For example, if the quenching processes in groups cause star 
formation to end on a time scale tq ~ 1 Gyr (whic h best 
reproduces the BO effect in IKodama fc Bowerl l200lT ). then 
DEEP2 groups in this picture would have begun quenching 
their member galaxies efficiently only 1 Gyr prior to the ear- 
liest epoch considered here, i.e., at a redshift z ~ 1.7 (for 
h = 0.7). 

This basic conclusion (that the quenching mechanism 
turns on only at z < 2) is consistent with observations 
of high star-formation rat es in massive galaxies at z ~ 
2 {e.g., (|Daddi et al.l 120041 )1. In a ddition, it is su p porte d 
by the stellar-population study of ISchiavon et alJ (2006), 
which shows that nonstar-forming galaxies in DEEP2 at 
z ~ 0.9 have mean stellar age s of roughly 1.5 Gyr (indica- 
tive of_coieirching at^z < 2) |Schiavon et al.l [20061 ). More- 
over lHarker et all (2006) have shown that a relatively simple 
model of stellar population evolution, in which star forma- 
tion is turned off starting at z ~ 2, can explain the build-up 
of the red sequence from z ~ 1 to z ~ 0. The conclusion that 
DEEP2 groups only started quenching their member galax- 
ies at z ~ 2 is also nicely consistent with the emerging theo- 
retica l picture of galaxy formation (e.g.. [Dekel fc Birnboiml 
2006), in which star formation is quenched efficiently only 
in haloes with masses above ~ 10 12 Mq at z < 2. We discuss 
this last point in more detail in § 17.31 below. 

It is reasonable to worry, though, that the contami- 
nation of the group population by interloper field galaxies 
might cause the two populations to appear indistinguish- 
able when, in truth, a tiny difference still exists. Certainly 
Figure [10] shows that such contamination weakens the ob- 
served difference between the two populations (although it 
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does not induce spurious evol utionary trends). H owever, the 
tests on mock catalogues in iGerke et all (|2005h show that 
the sample of DEEP2 galaxies in VDM groups is dominated 
by galaxies that are truly in groups, so it should be pos- 
sible to discern all but the tinie st differences betwee n the 
samples. In addition, recall that ICooper et alJ (|2006ah find 
a similar convergence in the red fraction of galaxies in the 
most underdense and most overdense environments. Since 
that study does not suffer from the same contamination ef- 
fects as the VDM group catalogue, we may be confident that 
the apparent convergence of the two populations' blue frac- 
tions cannot be ascribed solely to interloper contamination. 

7.3 Comparison to semi-analytic models of galaxy 
formation 

It will be instructive to briefly compare the basic picture 
that the observations in §[5] imply — namely, that groups are 
not efficient at quenching star formation at z ~ 2 — with the 
picture that emerges from semi-analytic models (SAMs) of 
galaxy formation. This comparison is especially worthwhile 
because most SAMs have been designed to reproduce the 
details of the galaxy population at z ~ 0, so comparisons 
to the high-z DEEP2 population will provide sharp tests 
of their correctness. We will limit ourselves here to qualita- 
tive descriptions and arguments, however, and defer detailed 
comparisons to future work. 

In making the comparison, we w ill sometimes use as 
a concrete example the model used bylCroton et al.l (|2006T l 
to populate the Millenium run (jSpringel et al.ll2005f ). a very 
large, high-resolution N-body simulation of cosmic structure 
formation. The elements of thi s model are qualitatively very 
similar to other recent SAMs l|Bowerll200i ICattaneo et al.l 
120061 ; iKang et al.ll2006l) . all of which have been successful 
at reproducing the observed luminosity functions and colour 
distributions of galaxies at z ~ 0. In the lCroton et al.l (|2006l ) 
model, dark matter haloes a re populated with two types of 
galaxy (|Springel et al.ll200ll ): each halo contains one central 
galaxy and may also incorporate satellite galaxies, whose 
parent haloes merged with the main halo at some time in 
the past (this conceptual division of the galaxy population is 
nearly universal in modern-day SAMs). Each halo is also as- 
sumed to contain a baryonic component consistent with the 
overall cosmic baryon fraction. This gas can fuel ongoing 
star formation if it cools and is accreted by the halo's cen - 
tral galaxy (e.q.. lWhite fc Reesll 19781 : 1 White fc Frenklll99lf ); 
hence, in order to cease star formation and permanently join 
the red sequence, a galaxy must stop accreting new cold gas 
onto its disk. 

Most SAMs admit two basic routes by which gas accre- 
tion can be stopped, associated with the two types of galaxy 
in the halo. In the case of satellite galaxies, their parent 
haloes' gas is simply added to the gas reservoir of the larger 
halo when they are accreted. Gas is then only allowed to 
accrete onto its central galaxy (which sits at the minimum 
of the potential well), so model satellite galaxies will cease 
star formation forever once they have exhausted the gas sup- 
ply in their discs. All satellite galaxies will thus eventually 
migrate to the red sequence via strangulation. For central 
galaxies, the situation is somewhat more complicated. In all 
haloes, infalling gas is shock-heated to the halo viral temper- 
ature, but, in low-mass haloes, that hot gas cools on very 



short timescales so that accretion onto the central galaxy 
is very rapid. Star formation never ceases permanently in 
such galaxies. In more massive haloes, the cooling time is 
longer, and the gas forms a quasi-static, hot halo of gas. 
The seperation of haloes into these two regimes — the "rapid 
coolin g" and "static hot halo" regimes — has a long history 
fe.o.. lRees fc Ostriker|[l977l ; iBlumenthal et alJll984h : recent 
simulations show that the regimes are divided by a rather 
sharp transition mass of ~ 3 x IO^Mq (|Birnboim fc Dekell 
1200a iKeres et alj|2004 ) . 

Regardless of the halo mass, however, hot gas near the 
centre can still cool, accrete onto the central object, and 
form stars, unless some other source of energy is available to 
stop the cooling; thi s is the well-known cooling-flow problem. 
ICroton et all (2006) propose to solve this problem via heat- 
ing from low-luminosity AGN lurking in the haloes' mas- 
sive central galaxies. They refer to this heating as "radio 
mode" AGN feedback, to distinguish it from winds created 
by bright AGN in the "quasar mode" (which should only 
temporarily eliminate cold gas from a galaxy). In essence, 
if the central supermassive black hole (SMBH) is massive 
enough, its average AGN energy output will be sufficient 
to offset cooling in the halo gas and stanch the flow of 
gas onto the central galaxy; t his basic idea has also been 
implemented in other SAMs dBowerl [2OO6I ; ICattaneo et al.l 
l200fj ; IKang et al.ll2006l). Since halo mass and SMBH mass 
are linked (e.g.. lHaring fc Rixll2004j ). this model effectively 
predicts a threshold mass above which accretion onto the 
central galaxy ceases; then, once the galaxy has exhausted 
its existing supply of cold gas, star formation ends. 

The threshold mass above which "radio mode" heat- 
ing can stop cooling flows is a few tim es 1O 12 M0 
i|Dekel fc Birnboim|[200g ; ICattaneo et al 1 120061 ; Croton et al. 
in prep.), roughly an order of magnitude smaller than the 
typical masses of DEEP2 groups. Thus, the Croton et al. 
picture implies that all galaxies in DEEP2 groups should 
eventually be quenched and join the red sequence. The only 
group members that will be blue (i.e., star-forming or re- 
cently star-forming) in that model will be those that recently 
fell in to their host group and have not yet exhausted the 
cold gas in their discs. If this picture is accurate, then the 
convergence of ft for the group and field populations shown 
in Figure [7] implies that a typical DEEP2 group member at 
z = 1.3 was not in a halo above the quenching mass until 
tq ~ 1 Gyr before that epoch. This is broadly consistent 
with hierarchical structure growth in a ACDM universe. In 
the Millenium Run, for example, one can trace the history of 
a typical DEEP2-group-like halo, with a mass of ~ 1O 13 M0 
at z — 1.3. On average, such a halo's most massive progen- 
itor at z ~ 2 had a mass of ~ 5 x 10 12 Mq — roughly the 
mass at which radio-mode heating becomes efficient. There 
arc 1.1ft -1 Gyr of cosmic time between z = 2 and z = 1.3, 
a reasonable timescale for quenching and migration to the 
red sequence. It is worth pointing out that this explana- 
tion is conceptually quite different than the BO effect: as 
we have used the phrase in this paper, the BO effect is sim- 
ply a consequence of the assumption that groups and clus- 
ters quench star formation (and thus it affects only infalling 
satellites) whereas the mechanism described here is the on- 



20 Brian F. Gerke et al. 



set of quenching in group-size haloes (which affects central 
galaxies as welljj. 

As we have seen, evolutionary results fit nicely, in a 
broad sense, with the current semi-analytic picture of galaxy 
formation. However, there are some difficulties in detail, es- 
pecially as regards satellite galaxies. Satellites in the model 
are unable to accrete new gas; such accretion is reserved 
for central galaxies, regardless of the parent-halo mass. Un- 
der these assumptions, groups at z ~ 1.3 should have a 
lower blue fraction than the field, unless all of their galax- 
ies were isolated at z > 2. Since this seems unlikely, it may 
be necessary to revise the assumption that all satellites are 
stripped of their halo gas with equal efficiency, regardless of 
host-halo mass. Indeed, there is alread y significant evidence 
that this is an oversimplification. The ICroton et al.l (120061 ) 
model significantly over-predicts the abundance of faint red 
galaxies (see their Figure 11). Also, there is evidence that 
a satellite galaxy's colour is correlated w ith the colour of 
its central galaxy (jWeinmann et af1l2006al ); this would not 
occur if all satellites wer e stripped of thei r gas w ith equal 
efficiency. Most recently, IWeinmann et all (|2006bl ) showed 
that the Croton et al. model severely underpredicts the blue 
fraction of SDSS satellites and argued that the efficiency 
with which gas is stripped from satellite galaxies must scale 
with the mass of the host halo. We note in passing that 
such a revision to the model would also account naturally 
for the correlation seen in Figure HJ] between ft and group 
richness, since more massive groups would quench star for- 
mation more efficiently. 

7.4 On the nature of the red field galaxies 

We close this section by commenting on the continued pres- 
ence of a significant red galaxy population in the field at 
high redshift, even when the group and field blue fractions 
have converged. One of our main conclusions in this work 
is that groups are the primary locus of quenching and red- 
galaxy formation, and the most naive interpretation of this 
statement would suggest that there should be no red galax- 
ies in the field. This is an oversimplification, of course, not 
only because isolated elli pticals are known to exist locally 
(e.g., IColbert et al.l l2001). but also since the red sequence 
will be conta minated by some fraction of du sty star-forming 
galaxies (e.q. lBell et al]|2004al ; [Weineril2005l '). One might also 
expect such galaxies to be more common at high redshift, 
when the average star-formation rate of the Universe was 
higher. It should therefore come as no surprise that DEEP2 
contains a non-negligible population of red field galaxies at 
all z, but it is not immediately clear whether dusty star- 
formers can account for all such galaxies. 

To address this issue, we make use of the HST /ACS 
imaging that has been obtained for a portion t he Extended 
Groth Strip subregion of the DEEP2 survey (|Davis et al.l 



8 Note, however, that this effect, on its own, will not be sufficient 
explain the growing number density of red galaxies over the entire 
DEEP2 redshift range. The evolution seen in Figure [7] appears 
to flatten at z < 1, so for the red sequence to continue its build- 
up, there must be an increase in the number density of groups. 
That is, the growth of the red galaxy population at z < 1 might 
still be attributable solely to the hierarchical growth of structure. 



2006). In this region there are eleven red field galaxies in 
the redshift range 1 .15 < z < 1.3. Five of these show spiral 
structure or other visible evidence of star-forming activity, 
while six appear to be nonstar-forming spheroidal or lentic- 
ular objects. Thus it appears that a significant population 
of early-type field galaxies remains at these redshifts; the 
red field population observed in DEEP2 cannot be fully ac- 
counted for by dusty star-formers. 

The issue is further complicated by group-finding er- 
rors. The VDM group-finder fails to identify some fraction 
of group members; these will appear in the field sample, 
and som e of them wil l be ea rly- types. However, the statistics 
given m iGerke et all (120051 ) imply that only 11% of the field 
galaxies identified in sample IV are actually group members 
(not the value of 6% quoted in § 12.21 since we reclassified a 
significant number of group galaxies as field galaxies in § [4j) . 
As shown in Table [TJ 12% of these— or 1.3% of the total 
field sample — will be red. Combining this with the observed 
fraction of dusty star-formers still does not account for the 
full population of red field galaxies (about 7% of the field in 
sample IV). 

It is therefore worth considering a few mechanisms by 
which bona fide early-type galaxies could enter the field 
population. Here the distinction between central and satel- 
lite galaxies will again be useful. The threshold mass for 
quenching in central galaxies, a few times 10 12 Mq, is some- 
what lower than the mass range of groups considered in this 
study; therefore, there should be a population of quenched 
central galaxies that do not have satellites that enter the 
DEEP2 sample; these would appear as red galaxies in the 
field. This picture naturally leads to the correlation between 
fb and limiting magnitude seen in Figure [5] The strength 
of this correlation appears to be the same both in groups 
and in the field. But if all galaxies above a given mass are 
quenched, and if halo mass correlates roughly with lumi- 
nosity, then brighter subsets of the galaxy population will 
naturally have lower blue fractions, both in groups and in 
the field. This explana t ion a lso fits nicely with the recent 
results of IConrov et al.l (2006), who show that the brightest 
galaxies' mass-to-light ratios have grown dramatically be- 
tween z = 1 and the present day while M/L has remained 
roughly constant for fainter galaxies — presumably because 
the brightest galaxies have ceased star formation but con- 
tinue to accrete new, non-luminous mass. 

In addition, it has been shown that winds from merger- 
driven quasar activity may expel cold gas fr om galactic disks 
|Springel et al.l 120051 : [Hopkins et aT] l2005h . If such a wind 
arises in a galaxy whose halo is less massive than the thresh- 
old at which "radio-mode" AGN feedback is important, then 
the galaxy should eventually reincorporate some of the ex- 
pelled gas and resume star formation. But if the timescale 
for reincorporation is longer than the lifetime of hot, blue 
stars (~ 1 Gyr), then such galaxies will temporarily join 
the red sequence. Such temporarily red galaxies could also 
account for some red field objects. 

It is also noteworthy that the blue fraction in the field 
appears to evolve at z > 1, if less strongly than in groups. 
This can be seen in the lower right panel of Figure [7] fb in 
the field changes by ~ 0.05 between z = 1.0 and z — 1.3, 
while the group value changes by ~ 0.1. It is worth explor- 
ing whether this effect is real, or whether it can simply be 
attributed to contamination of the field sample by group 
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galaxies, arising from group-finding errors. Let us assume, 
to start, that the true blue fraction in the field is constant 
with z. Then, if the field sample contains a fraction Cf of 
contaminating group galaxies, and the group sample simi- 
larly has a contamination fraction c 9 of field galaxies, it is 
straightforward to show that, for a given observed change 
A/^ roup in the group blue fraction, the observed change in 
the field blue fraction is 

A/ flcld = ^_ A/ group^ (1Q) 

1 Cg 

For sa mple IV, the success statistics given in iGerke et al.l 
|2005l ) imply that c g — 0.46 and c/ = 0.11 (as discussed ear- 
lier in this Section). Thus, the observed change in the field 
blue fraction should be ~ 20% of the observed change in the 
group blue fraction. This ratio is, in fact, marginally con- 
sistent with the error bars in Figure [7] Hence, the observed 
decline in the field blue fraction with time is largely caused 
by contamination by misclassified grou p galaxies. This con - 
clusion is supported by the findings of lCooper et al.l (|2006al ). 
who observe no evolution in the red fraction of galaxies in the 
most underdense environments. On the basis of the present 
study alone, however, it remains possible that there is also 
some small amount of intrinsic evolution in the field popu- 
lation. Such evolution would not be too surprising: since the 
characteristic halo mass M* does not reach th e quenching 
mass of a few times 10 12 A/q until z ~ 1 (e.g., Ree d et ah! 
2003), the number density of quenched central galaxies in 
the field will be growing relatively rapidly at earlier times. 
Regardless, it remains clear from Figures [7] and [9] that the 
build-up of the red sequence in DEEP2 takes place prefer- 
entially in groups. 



8 CONCLUSION 

We ha ve used the DEEP2 group catalogue of IGerke et al.l 
(2005) to study the blue fraction of galaxies in groups and 
the field at z ~ 1. After creating four samples that are care- 
fully designed to be free from colour-dependent and redshift- 
dependent selection effects, we probe the dependence of ft 
on a galaxy's environment (group or field), on galaxy lumi- 
nosity, on group properties, and on redshift. 

First, over the redshift range 0.75 z < 1, there is a sig- 
nificant difference in ft between the DEEP2 group and field 
populations, with a substantially lower in groups than in 
the field. The well-known colour segregation observed be- 
tween nearby group and field galaxy populations is thus al- 
ready in place at z ~ 1. In addition, there is a negative 
correlation between /;, and group richness over this redshift 
range, echoing local correlations between galaxy colour and 
various proxies for halo mass. A trend also exists at these 
redshifts between and the luminosity of the galaxies be- 
ing considered, with fainter galaxies being more likely to be 
blue; the strength of this trend does not appear to depend 
on whether those galaxies are in groups or the field. 

No siginificant evolution in ft occurs in the DEEP2 
group or field galaxy populations over the range 0.75 z < 
1; however, at z > 1, there is dramatic, statistically signifi- 
cant evolution in the group population, with fi rising rapidly 
as z increases. Most notably, the value of ft in groups be- 
comes indistinguishable from the field value at z ~ 1.3. To 
assess the robustness of these results, we have constructed 



mock catalogues that replicate the colour-environment cor- 
relations in the DEEP2 data. Tests with these mock cata- 
logues indicate that the evolution and convergence are not 
artifacts of the group-finding algorithm. Moreover, these 
tests suggest that the observed evolution can be explained 
partially — but not entirely — by evolution in the distribution 
of galaxy environments, i.e. by the growth of large-scale 
structure, coupled with the known correlation between en- 
vironment and galaxy properties. However, to explain the 
strong evolution and convergence, some further mechanism 
is required. 

Because the ft evolution appears to weaken at z < 1, 
we conclude that it is different in nature from the classical 
BO effect, in which fb in clusters rises dramatically from 
z — out to z ~ 0.5. Indeed, the BO effect is a natural 
outcome of hierarchical structure formation, provided that 
groups and clusters are already efficient at quenching star 
formation in the galaxies that fall into them. But there is a 
convergence between ft, values in groups and in the field at 
z ~ 1.3, which implies that DEEP2 groups were, at some 
epoch, not especially efficient at quenching infalling galaxies. 
In particular, if it takes ~ 1 Gyr for a quenched galaxy to 
join the red sequence, our results indicate that the groups in 
this study only became suitable environments for quenching 
at z ~ 2. 

This conclusion is broadly consisten t with current semi- 
analyt i c mo d els of galaxy formati o n (ICroton et alj [20061 ; 
iBoweij l200d : ICattaneo etafl 120061 ; iKang et al.1 120061 ). In 
those models, the central galaxies of dark matter haloes 
are quenched when low-level AGN feedback becomes strong 
enough to stop gas in the halo from cooling and condens- 
ing onto the galaxy to form stars; this process becomes 
efficient at a threshold halo mass of a few ti mes 1O 12 M0 
(|Dekel fc Birnboimll2006l ; ICattaneo et~ai1l2006l ; Croton et al 
in prep.). Our mock catalogues imply that the groups ob- 
served in DEEP2 have masses in the range 5 x 10 12 Mq < 
M < 5 x 1O 13 M0. A typical halo in this mass range at 
z ~ 1.3 will drop below the threshold mass for quenching 
~ l/i _1 Gyr earlier, at z ~ 2. This is a potential explana- 
tion for the convergence in the group and field ft values at 
2 ~ 1.3. Present-day SAMs may have difficulty explaining 
these observations in detail, however, particularly because 
they universally assume that satellite galaxies are quenched 
in all haloes regardless of mass or redshift. In the future, it 
will be worthwhile to make detailed comparisons between 
our observations and the precise predictions of SAMs, in or- 
der to refine the physical description of galaxy formation 
that they provide. 
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